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DESIGN  OF  SOLAR  HEATING 
AND  COOLING  SYSTEMS 


1 INTRODUCTION 


Background 

The  technical  feasibility  of  heating  and  cooling 
buildings  using  Hat-plate  solar  collectors  has  been 
established  both  in  theory  and  practice.  Although  im- 
provements in  component  and  system  design  are  forth- 
coming. one  can  approach  solar  heating  and  cooling 
technology  with  full  confidence  that  a practical,  reliable 
system  can  be  constructed.  The  construction  phase  of 
solar  systems  requires  little  more  skill  than  is  presently 
required  to  install  conventional  heating  and  cooling 
systems;  the  design  phase,  however,  is  considerably 
more  complex. 

Sunlight  is  inherently  an  intermittent  source  of 
energy  at  the  earth’s  surface.  For  this  reason,  solar 
energy  systems  for  heating  domestic  hot  water  or  for 
heating  and  cooling  buildings  are  not  and  frequently 
cannot  be  designed  to  meet  the  full  demands  of  the 
building  or  buildings  being  served.  System  design  is 
therefore  uncoventional  in  that  solar  energy  supplies 
some  but  not  all  of  the  building’s  energy  requirements. 
The  fundamental  problem  in  analyzing  solar  energy 
systems  is  determination  of  the  collector  array  area 
which  will  provide  the  greatest  cost  benefit  in  meeting 
the  highest  fraction  of  the  building’s  energy  load.  In 
the  past,  because  of  the  uncontrolled  nature  of  solar 
energy,  extensive  and  costly  computer  studies  (taking 
into  account  hourly  weather  data  from  the  site  in 
question)  have  been  required  to  evaluate  proposed 
solar  projects  and  to  design  solar  energy  systems.  The 
tools  used  for  these  studies  have  generally  been  propri- 
etary and  thus  not  available  to  Corps  of  Engineers 
Districts.  Those  which  have  been  available  have  generally 
been  difficult  to  use  and  are  not  ready  for  general  use 
by  design  engineers. 

Hence,  there  is  a basic  need  for  a simple,  convenient 
method  for  making  the  necessary  design  calculations 
and  for  assessing  energy  and  economic  benefits.  In 
particular,  simplified  techniques  for  making  preliminary 
energy  and  economic  analyses  are  needed  to  minimize 
the  number  of  expensive  computer  studies  required. 
Where  computer  studies  are  required,  user-oriented 
computer  programs  are  necessary. 


Objective 

The  objective  of  this  study  was  to  develop  and  field 
test  a straightforward  procedure  that  would  allow 
District  Engineers  to  perform  energy  and  cost  benefit 
analysis  as  well  as  final  design  calculations  for  individual 
solar  energy  system  applications.  The  procedure  would 
permit  ( I ) prelimmuiy  solar  energy  system  performance 
assessment  from  charts  and  graphs,  and  (2)  the 
computer-aided  evaluations  of  solar  energy  system  per- 
formance necessary  for  final  design.  The  procedure  was 
to  be  applicable  to  solar  domestic  hot  watei.  solar 
heating  only,  and  to  heating  and  cooling  systems.  This 
report  describes  how  to  use  the  graphical  and  computer 
tools  developed  in  response  to  the  study  objective. 

Approach 

The  approach  used  in  this  research  consisted  of  the 
following  steps: 

1.  Develop  a computer  simulation  progiam  lor 
analyzing  the  performance  of  solar  energy  systems. 

2.  Use  the  solar  simulation  program  to  determine 
the  performance  of  typical  solar  systems  in  meeting 
predetermined  hourly  building  loads  at  each  of  five 
sites  in  the  continental  United  States  using  hourly 
incident  radiation  data. 

3.  Study  the  effects  on  solar  system  performance 
of  varying  solar  system  design  parameters  (including 
collector  area,  collector  tilt  angle,  storage  tank  volume, 
and  heat  exchanger  effectiveness). 

4.  Based  on  the  results  of  the  parametric  studies, 
develop  a graphical  method  for  estimating  solar  system 
performance.1-2 

5.  Develop  a user-oriented  computer  simulation 
program  for  making  rigorous  solar  system  performance 
calculations  during  the  final  design  phases.3 

'll.  C.  Hittle.  D.  I-.  Holshouser.  and  C>.  N.  Walton. Interim 
Feasibility  Assessment  Guidance  for  Solar  Heating  and  Cooling 
of  Army  Buildings.  Technical  Report  E-91/ADA026588  (l  .S. 
Army  Construction  Engineering  Research  Laboratory  |CT.RL|, 
May  1976). 

2D.  C.  Hittle.  G.  N.  Walton,  0.  I . Holshouser,  and  1).  J. 
Leverenz,  Predicting  the  Performance  of  Solar  Energy  Systems. 
Technical  Report  E-98/ADA035608  (CTRL.  January  1977). 

3D.  C.  Hittle.  Building  Loads  Analysis  and  System  Thermo 
dynamics  (BI.ASTI  Program  Vo I I.  Lscrs  Manual.  Technical 
Report  E-U9/ADA048734  (GI  RL.  November  1977):  and 
Building  Loads  Analysis  and  System  Thermodynamics  (BI.A  ST) 
Program:  Vol  II.  Reference  Manual.  Technical  Report  1-1  19/ 
ADA048982  (CF.RL,  December  1977). 
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o Develop  lecoutmended  procedures  lor  using  lire 
graphical  and  computer  simulation  ieelini<|ues  lor 
making  solar  system  energy  and  economic  analyses, 
and  for  optimizing  solar  system  designs  based  oil  the 
highest  benetit-to-cost  ratio. 

7.  Field  lest  the  graphical  and  computer  procedures 
for  making  a solar  system  energy  and  economic  analysis 
at  a Corps  District  to  insure  that  the  procedures  are 
eleai  and  user-oriented,  and  that  the  original  input  data 
are  available. 

X.  Publish  the  graphical  and  computer  procedures, 
modified  after  the  field  test,  as  an  engineer  Technical 
Letter. 

Organization  of  Report 

Chapter  2 introduces  the  graphical  method  for  com- 
puting expected  performance  of  solar  collector  arrays, 
and  Chapter  3 presents  examples  of  its  use.  Chapter  4 
describes  a method  for  optimizing  the  solar  energy 
system  on  the  basis  of  highest  benefit-to-cost  ratio. 
Chapter  5 describes  the  computer-aided  method  for 
solar  energy  systems  which  is  recommended  for  per- 
forming design  calculations  on  larger  projects.  Chapter 
6 outlines  some  practical  considerations  relevant  to  the 
application  of  solar  technology. 

Mode  of  Technology  Transfer 

The  evaluation  procedures  and  information  in 
Chapters  2 through  6 will  be  furnished  to  District  and 
Facilities  engineers  as  an  engineer  Technical  Letter. 


THE  UNIVERSAL  SOLAR  SYSTEM 
2 PERFORMANCE CURVE-A 
GRAPHICAL  APPROACH 


Introduction 

During  the  development  of  a solar  system  simula- 
tion model  and  design  methodology,  the  U.S.  Army 
Construction  Engineering  Research  Laboratory  (CERL) 
performed  several  hundred  solar  system  simulations  for 
typical  Army  buildings  in  various  paris  of  the  country. 
Analysis  of  the  solar  system  performance  curves  for 
these  systems  indicated  that  with  proper  normalization, 
the  performance  of  a given  solar  system  for  the  various 
buildings  in  all  locations  could  be  represented  by  a 
single  universal  performance  curve. 


I igme  I is  a schematic  ol  the  type  <d  system  studied 
in  developing  llie  universal  cuive.  SoIji  ladialion.  when 
available,  is  converted  to  thermal  energv  at  the  col- 
lector and  storage  pumps.  The  heat  exchanger,  which 
isolates  the  collection  and  storage  loops,  allows  the 
collector  iluid  to  be  freeze-protected.  Normal  operation 
of  this  system  permits  heating  of  the  storage  tank 
whenever  sufficient  solar  energy  is  available. 

Heal  is  taken  from  the  tank  by  the  load  pump 
whenever  there  is  a demand  foi  heat,  whether  it  is  at  a 
healing  coil  or  al  the  geneiatoi  of  an  absorption  an 
conditioner.  As  shown,  this  system  is  Representative  ol 
a large  class  of  liquid  solar  energy  systems,  and  all  the 
simulations  run  al  CERL  weie  based  on  such  a system. 

The  universal  curve  may  be  consulted  to  determine 
the  overall  performance  of  the  system  pictured  in 
Figure  1 . Three  such  curves  are  presented  in  litis  report : 
one  for  solar  domestic  hot  water,  one  lor  solar  heating 
only,  and  one  for  solar  healing  and  cooling  applications. 
It  should  bo  pointed  out  that  the  development  ol  these 
curves  was  based  on  several  assumptions,  l irsl.  a 
single-cover,  fiat-plate  collector  with  an  emissivity  of 
. 10  and  an  absorptivity  of  .00  was  used  as  tbit'  reference 
collector.  The  tilt  of  these  collector.,  fro  -,  the  hori- 
zontal depended  on  the  particular  application.  For 
solar  domestic  hoi  water,  this  angle  was  assumed  to  be 
equal  to  the  latitude;  for  solar  space  heating,  the  angle 
was  taken  as  the  latitude  plus  10  degrees;  for  solar 
heating  and  cooling,  it  was  latitude  minus  10  degrees. 
These  collection  tilt  angles  were  found  to  be  near 
optimum  lor  these  applications.  All  collectors  were 
unshaded  and  facing  south.  The  heal  exchanger  ef- 
fectiveness was  laken  to  lie  .X,  and  the  storage  tank 
held  16  lb  of  water  per  square  foot  of  collector  (78 
kg/m2  or  about  2 gal/sq  ft).  The  tank  insulation  had  an 
“R"  value  of  |0.  (These  values  of  heat  exchanger 
effectiveness,  tank  volumes,  and  lank  insulation  were 
considered  reasonable  for  this  type  of  system.)  Line 
losses  were  assumed  to  be  negligible  in  all  cases 

While  variations  in  all  these  parameters  will  affect 
system  performance  (see  Appendix  A),  ii  was  found 
that  the  assumed  collector  efficiency  is  by  far  the 
dominant  (actor.  A method  for  adjusting  collector  area 
for  types  other  than  the  reference  collector  is  presented 
in  Chapter  3.  Variations  in  the  other  parameters  can  be 
evaluated  in  the  final  design  phase  using  the  computer 
simulation  model  (see  Chapter  5).  Since  one  purpose  of 
this  report  is  to  present  a simple  method  to  be  used  for 
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Figure  1.  Schematic  of  the  liquid  solar  energy  system  used 
in  the  development  of  the  universal  curves. 


feasibility  studies,  explicit  correction  factors  for  the 
effects  of  varying  parameters  other  than  the  collector 
type  have  been  excluded. 

Use  of  the  universal  curves  requires  calculation  of 
building  energy  loads  and  compilation  of  on-site 
weather  data.  Once  this  information  is  available,  the 
curve  relates  solar  system  performance  to  the  collector 
array  area.  In  order  to  apply  the  curve,  however,  an 
understanding  of  the  following  terms  is  required: 

I.  Annual  or  monthly  incident  solar  radiation  flux 
density.  0C  This  is  the  solai  flux  density  on  the  tilted 
collector  array  in  Btu/square  foot/month  or  Btu/square 
foot/year  (langleys/month  or  langleys/year).  Solar 
radiation  data  are  normally  given  in  terms  of  horizontal 
radiation  densities.  Thus,  these  numbers  must  be  cor- 
rected for  the  tilt  angle  of  the  collector  array.  The 
method  for  making  these  corrections  is  presented  in 
the  Solar  Radiation  Data  section  of  this  chapter. 


[Hq  1 1 


where  Ac  is  the  collector  area.  This  quantity  is  essen- 
tially a measure  of  the  available  solar  energy  (incident 
on  the  collector  array  of  the  system  in  question)  com- 
pared to  the  total  energy  required. 


4.  Percentage  of  the  energy  requirements  met  by 
solar,  p.  This  term  is  the  fraction  of  QL  met  by  the 
solar  energy  system.  The  quantity  p may  be  calculated 
from  the  appropriate  universal  curve  of  solar  system 
performance.  Figure  2 shows  the  form  of  this  curve  for 
solar  domestic  hot  water.  From  the  figure,  it  can  be 
seen  that  the  relationship  between  Ps  andp  allows  cal- 
culation of  the  solar  system  performance.  For  domestic 
hot  water. 


2.  Annual  and  monthly  energy  requirements,  Q( . 
For  the  purposes  of  the  universal  curve,  these  energy 
requirements  are  defined  as  the  total  thermal  energy 
in  the  form  of  hot  water  required  by  the  building’s 
domestic  hot  water  heating  system,  space  heating 
system,  or  absorption  chiller. 

3.  Solar  system  performance  parameter,  P . The 
solar  system  performance  parameter  is  defined  as  the 
ratio  of  annual  (or  monthly)  radiation  incident  on  the 
collector  array  to  the  annual  (or  monthly)  energy  re- 
quirements of  the  building,  as  given  by  Eq  1 : 


Closed  form  expressions  of  the  universal  curve  for  solar 
heating  and  cooling  and  solar  heating  only  are  given  in 
Chapter  3. 

As  an  example  of  the  application  of  the  method, 
consider  a 200-man  barracks  in  Topeka.  KS.  which 
requires  1.73  X 109  Btu  (1.83  X 109  kJ)  of  energy 
annually  to  heat  hot  water.  Use  of  the  universal  curve 
permits  determination  of  the  percentage  of  this  load 
which  can  be  met  by  2400  sq  ft  (223  m2)  of  collectors 
tilted  at  39  degrees  from  the  horizontal. 


annual  incident  radiation /annual  energy  requirements,  p* 


The  solar  radiation  llux  density,  0 , on  a collector 
tilted  at  39  degrees  is  6.0  X 105  Btu/sq  ft  (6.8  X I06 
kj/m2).  (The  method  for  obtaining  solar  flux  density 
on  a tilted  surface  is  described  in  the  Solar  Radiation 
Data  section.)  From  Kq  1 , P = .83;  from  Eq  2,  p is  .48. 
Thus.  2400  sq  ft  (223  in2)  of  collector  would  meet  48 
percent  of  the  load;  i.e.,  the  collector  array  would  sup- 
ply (.48X1.73  X 109  Btu)  = 8.3  X 108  Btu(8.8X  108 
kj),  and  the  auxiliary  system  wouid  have  to  supply  the 
remaining  9.0  X 108  Btu  (9.5  X 108  kJ)  of  hot  water 
by  means  of  a conventional  hot  water  heating  system. 

If  the  efficiency  of  the  auxiliary  energy  supply  sys- 
tem (an  oil  boiler,  for  example)  is  80  percent,  the  solar 
energy  system  would  save  1.04  X I09  Btu  ( 1 . 1 0 X I0!> 
kJ ) of  fuel  energy.  If  the  heating  value  of  oil  is  1 38.000 
Btu/gal  (38.464  kJ/l !)  at  a cost  of  $0.50/gal  (SO.14/12). 
the  annual  dollar  savings  would  be  $3,768.  This  savings 
could  be  used  to  offset  the  capital  cost  of  the  solar  sys- 
tem. With  this  information,  an  economic  cost  benefit 
analysis  of  this  candidate  solar  system  can  be  made. 
The  procedure  can  be  repeated  for  different  collector 
areas  until  the  collector  area  giving  the  maximum  cost 
benefit  is  found. 

It  should  be  pointed  out  that  there  is  another  class 
of  solar  systems  which  operate  by  wanning  and  circulat- 
ing air  through  the  collectors  rather  than  water.  These, 
however,  are  not  as  widely  used  as  liquid  solar  energy 
systems.  (Air  systems  have  been  used  with  some  success 
in  residential-scale  solar  domestic  hot  water  and  heating- 
only  applications,  but  are  not  suited  for  solar  air- 
conditioning.)  A study  comparing  the  liquid  and  air 
systems  has  shown  that  for  these  heating  applications, 
the  performance  of  the  two  is  almost  identical.4 

Obtaining  Required  Input  Information 

As  the  example  in  the  previous  section  shows,  two 
sets  of  input  information  in  addition  to  collector  area 
are  required  to  use  the  universal  curve;  the  thermal 
energy  demand.  Q,,  and  the  radiation  tlux  density, 
Oc.  The  following  sections  describe  how  this  infortna- 
iion  can  be  obtained. 

liiergy  Consumption  Data 

The  type  of  data  required  to  determine  QL  and  the 
sources  of  that  data  depend  on  the  type  of  system  being 
considered.  Since  solar  energy  systems  produce  thermal 


xieorge  O.  G.  LSI',  “Heating  of  Buildings  with  Solar 
Energy."  Solar  Heating  and  Cooling  of  Buddings,  Vol  3,  Joint 
Conference  Proceedings  of  American  Section,  International 
Solar  Energy  Society  and  the  Solar  Energy  Society  of  Canada, 
K.  \V.  Boer,  ed.  (August  15-20,  1976). 


energy  in  the  form  of  hot  water,  consumption  or  energy 
requirements  in  any  other  terms  such  as  kilowatt  hours 
supplied  to  or  required  by  a centrifugal  chillei.  Bin's 
of  chilled  water  required,  gallons  of  oil  or  cubic  feet 
of  gas  supplied  to  a boiler,  or  gallons  of  hot  water 
used-must  be  converted  to  thermal  energy  in  the  form 
of  hot  water  (or  hot  air  if  a hot  air  system  is  used.) 
This  usually  requires  taking  into  account  the  efficiency 
of  conversion  equipment  as  well  as  the  appropriate 
conversion  factors.  Determination  of  for  the  three 
types  of  systems  is  described  below. 

Solar  Domestic  Hot  Water  System.  To  use  the  uni- 
versal curve  for  domestic  hot  water  systems,  the  annual 
energy  needed  for  heating  the  hot  water  is  required. 
For  the  purpose  of  designing  a solar  energy  system,  the 
best  source  of  this  information  is  metered  data  from 
the  building  in  question  or  from  an  identical  or  similar 
building  at  the  same  location.  Metered  consumption 
data  for  either  energy  or  fuel  input  to  the  hot  water 
heater  or  gallons  of  hot  water  consumed  can  be  used. 
When  fuel  consumption  is  known,  corrections  for  boiler 
efficiency  must  be  taken  into  account  as  discussed 
above.  When  actual  consumption  data  cannot  be  ob- 
tained. data  can  be  derived  from  various  handbooks 
(such  as  the  section  “Plumbing”  in  TM  5-810-1 5 ) or 
from  local  or  national  plumbing  codes.  However,  care 
must  be  exercised  in  using  data  from  these  sources, 
since  they  are  usually  developed  for  design  purposes 
and  thus  are  often  more  closely  related  to  peak  demand 
than  to  average  hot  water  consumption.  Whichever 
method  is  used,  if  annual  hot  water  demand  is  deter- 
mined, the  annual  energy  consumption  is  obtained 
from  Eq  3: 

Q,  = DwC  (T  -T.  ) |Eq  3] 

L p'  out  in'  1 “ I 

where  QL  = annual  energy  consumption  in  Btu(J) 

D = annual  hot  water  demand  in  gal  (V) 
w = density  of  water  = 8.33  lb/gal  ( 1 .00  kg/?) 

C = specific  heat  of  water  = 1 Btu/lbm  °F 
P (4.1868  kJ/kg°K) 

T = temperature  of  hot  water  supply  in  °F 
<°C) 

Tjn  = temperature  of  supply  water  in  °F  (°C) 

Solar  Heating  Systems.  Solar  energy  incident  upon  a 
collector  array  in  the  summer  is  of  little  use  for  solar 
space  heating  the  following  winter.  For  this  reason,  the 


f Mechanical  b gn:  Heating.  Ventilating  and  Air  Condi 
Honing  Changes  1-5.  TM  5-810-1  (Department  ot  the  Army. 
January  1956). 
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universal  curve  for  solar  heating  must  be  applied  on  a 
monthly  basis  only  for  those  months  when  a space 
heating  load  occurs.  This,  in  turn,  requires  that  the 
monthly  energy  demand  for  the  healing  be  obtained. 

Once  again,  the  best  source  of  data  for  estimating 
the  amount  of  heat  energy  required  for  a particular 
building  or  set  of  buildings  is  measured  data  for  the 
building  in  question  or  for  similar  or  identical  buildings 
at  the  same  location.  These  data  may  only  be  available 
m terms  of  the  amount  of  fuel  oil  used  per  month  or 
the  amount  of  gas  consumed  in  heating  the  particulai 
building  or  buildings.  It  so,  the  fuel  consumption  data 
must  be  adjusted  to  account  for  the  efficiency  ol  the 
boiler  01  furnace  supplying  the  heat  in  the  space.  Typi 
callv.  solar  heating  systems  include  both  heating  the 
building  and  the  domestic  hot  water  for  the  building. 
This  can  be  accounted  for  in  the  energy  requirement 
by  adding  the  monthly  domestic  hot  water  load  to  the 
heating  load. 

If  measured  data  for  the  building  being  considered 
are  not  available,  the  building's  heating  energy  demand 
can  be  estimated  using  one  of  several  methods.  One  such 
method,  described  in  a CURL  report  entitled  Energy 
Utilization  Index  Method  for  Predicting  Building 
Energy  Use ,6  provides  a manual  procedure  to  compute 
monthly  and  annual  heating  and  cooling  loads.  This 
procedure  was  developed  from  extensive  computer 
simulations  and  uses  readily  available  building  descrip- 
tions and  weather  data  as  its  input. 

The  Carrier  Rational  Hnergy  Analysis  Procedure 
I RLAP)7  is  another  method  for  calculating  building 
heating  requirements.  This  is  a modified  bin  method8 
which  requires  the  calculation  ofheating  loads  at  several 
different  design  points  of  the  system.  These  loads, 
combined  with  frequency  of  occurrence  of  weather 
conditions,  which  can  be  obtained  from  Air  Force 
Manual  88-8,9  provide  an  estimate  of  the  monthly 
energy  requirements  for  heating. 


6L.  M.  Wuwlinidiind  and  1).  C.  Hittle,  Hnergy  Utilization 
Index  Method  for  Predicting  Building  Energy  Use,  Vols  I & 11 
Interim  Report  F-105/ADA039913  and  ADA040344  (CURL. 
May  1977). 

1 Rational  Energy  Analysis  Procedure  (REAP)  (Carrier  Air 
Conditioning  Company,  undated). 

*ASHRAI:  Handbook  of  Eundamentals  (American  Society 
of  Heating.  Refrigerating  and  Air  Conditioning  Engineers 
IASHRAI  |,  1972). 

9 Engineering  Weather  Data.  Al'M  88-8,  Chapter  6 (Depart- 

ment of  the  Air  Force.  1976). 


Another  source  of  loads  data  is  the  output  from 
energy  analysis  computer  programs  such  as  the  Building 
Loads  Analysis  and  System  Thermodynamics  ( BLAST ) 
program.10  Loads  information  obtained  from  a com- 
pute! simulation  of  the  building  under  construction 
could  be  used  for  the  solar  system  analysis.  However,  it 
such  a simulation  had  not  been  performed  during  the 
design  phase,  the  additional  effort  required  to  prepare 
a building  input  deck  would  probably  not  be  warranted 
for  a solar  energy  feasibility  study  alone 

Montly  energy  consumption  foi  domestic  hot  water 
heating,  which  is  generally  included  in  a solai  heating 
system,  can  be  determined  by  applying  the  methods 
described  in  the  preceding  section  on  a montly  basis.  It 
should  be  noted  that  measured  data  often  have  domestic 
hot  water  and  space  heating  combined. 

Design  day  calculations  alone  cannot  be  used  to 
obtain  heating  loads,  since  they  provide  peak  energy 
demands  rather  than  (lie  average  demand  requned  for 
estimating  monthly  energy  consumption.  Therefore, 
the  methods  for  determining  heating  loads  described  in 
the  ASIIRAL  Handbook  oj  Eundamentals  or  in  the  De- 
partment of  the  Army's  TM  5-810-1  should  be  avoided. 

Solar  Heating  and  Cooling  Systems.  The  universal 
curve  for  solai  heating  and  cooling  systems  is  an  annual 
curve,  since  there  is  a year-round  demand  for  energy. 
The  Ql  required  is  the  sum  of  the  heating  and  cooling 
energy  requirements  (plus  the  domestic  hot  water  heat- 
ing energy  requirement  if  it  is  to  be  included  in  the 
system). 

The  annual  energy  required  for  heating  and  cooling 
is  again  best  determined  from  measured  data  for  the 
building  in  question  or  a similar  or  identical  building  at 
the  same  location.  In  determining  the  energy  required 
for  cooling,  attention  must  be  given  to  the  differences 
in  the  coefficients  of  performance  (COP)  of  the  vaiious 
chillers  involved.  For  example,  if  measuring  the  con- 
sumption of  electricity  for  a conventional  chiller  in  the 
building  indicates  that  1.000.000  kWh  are  used  annual- 
ly, then  this  number  must  first  be  converted  to  Bin's 
by  multiplying  the  8412,  giving  3.41  X 10s  Btu  (or 
converted  to  kilojoules  by  multiplying  by  3o00.  giving 
3.60  X 109  kJ),  and  then  multiplied  by  th  e centrifugal 
chiller’s  mean  COP.  If  the  mean  COP  is  4.  13.6  X 109 
Btu/yr  ( 14.3  X I09  kJ/yr)  are  required  to  be  delivered 


l0D.  C.  Hittle.  Buihtwg  t.oads  Analysis  and  System  thermo 
dynamics  (BI.ASTl  Program  Users  Manual.  Technical  Report 
E-1 19/ADA048734  (CURL,  December  1977). 
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by  the  chiller  to  the  cooling  system.  This  chilled  water 
demand,  when  divided  by  the  COP  of  a solar  energy 
system’s  absorption  chiller  ( say  .05)  indicates  that 
21.0  X 109  Btu  (22.2  X I09  kJ)  must  be  supplied  to 
the  absorption  chiller  to  meet  the  annual  energy  re- 
i|uired  for  cooling.  Since  heating  and  cooling  systems 
may  also  include  domestic  hot  water  heating,  total 
demand  for  heat  energy  for  heating  domestic  hot  water 
can  be  summed  with  the  total  demand  for  heating  and 
added  to  the  annual  demand  for  cooling  energy  in  order 
to  determine  the  annual  building  energy  requirement, 

Ut- 
il' measured  consumption  data  are  not  available,  the 
energy  use  index  method,  REAP  method,  or  BLAST 
simulation  program  described  in  the  previous  section 
can  be  used  to  determine  annual  heating  and  cooling 
loads.  Again,  calculations  based  on  peak  cooling  loads 
should  not  be  performed. 

Solar  Radiation  Data 

The  tilt  angle  of  the  collector  plate  greatly  affects 
the  amount  of  solar  radiation  striking  the  surface  of 
the  plate.  Most  measured  and  reported  data  are  for 
solar  radiation  incident  upon  a horizontal  surface.  This 
section  describes  where  to  obtain  this  information  and 
how  to  correct  it  for  tilted  collector  arrays. 

Average  solar  radiation  data  for  many  sites  around 
the  country  are  published  in  the  National  Oceanic  and 
Atmospheric  Administration's  (NOAA)  Climatic  Atlas 
of  the  United  States.  Annual  and  monthly  solar  radia- 
tion maps  from  this  document  are  reproduced  in 
Appendix  B.  Since  the  radiation  values  on  these  maps 
are  given  in  daily  means,  the  annual  and  monthly  values 
are  obtained  by  multiplying  by  the  number  of  days  in 
the  year  or  month,  depending  on  the  map  being  used. 

The  Climatic  Atlas  of  the  United  States  also  contains 
tabulated  summaries  of  radiation  data  for  specific  sites. 
More  detailed  summaries  are  frequently  available  from 
NOAA  or  from  local  weather  services  or  state  agencies. 
Data  for  the  particular  site  in  question  should  be  ob- 
tained, if  possible,  since  incident  solar  radiation  fre- 
quently varies  substantially  over  relatively  small  geo- 
graphical distances.  This  is  particularly  true  in  coastal 
regions  and  in  regions  near  mountains,  where  local 
climatological  variations  are  severe. 

Once  the  collector  tilt  angle  is  known,  the  radiation 
flux  density  on  the  tilted  collector  surface  can  be  deter- 
mined from  horizontal  radiation  data.  Fortunately, 
there  is  an  optimum  collector  tilt  angle  (measured  from 


horizontal)  for  each  type  of  system.  For  heating  and 
cooling  systems,  this  optimum  tilt  angle.  i>r,  is  roughly 
the  location  latitude  minus  10  degrees,  foi  heating-only 
systems,  it  is  roughly  the  latitude  plus  10  degrees:  and 
for  domestic  hot  water  systems,  it  is  roughly  equal  to 
the  latitude.  These  angles  provide  lot  collection  ol  the 
greatest  amount  of  solar  energy  in  each  application. 
Variations  of  t 5 degrees  affect  performance  only  sligli  i- 
ly.  In  all  cases  the  optimum  azimuth  angle  is  due  south 
and  again,  slight  deviations  from  due  south  (±  10  de- 
grees) do  not  significantly  reduce  system  performance. 

Once  the  optimum  tilt  angle  lias  been  determined, 
the  annual  or  monthly  radiation  flux  density.  Qc<  on 
the  optimally  tilted  surface  can  be  estimated  from  Eq4: 

QC  = KH„  |Eq4] 

where  K = a correction  factor  depending  on  the  col- 
lector tilt 

Fly  = the  annual  or  monthly  radiation  flux 

density  on  a horizontal  surlace  (Btu/sq  ft 
or  J/m2). 

For  solar  systems  for  which  annual  solar  radiation 
data  are  used  (i.e.,  heating  and  cooling  or  domestic  hot 
water),  the  correction  factor,  K,  is  given  approximately 
by 

cos(u  -7  -v  ) 

K = -^-7)  ^ 

where  vL  = the  latitude  in  degrees 

uc  = the  optimum  angle  from  the  horizontal 
in  degrees. 

Eq  5 is  an  empirically  derived  equation  based  on  th<* 
results  of  simulation  studies  and  is  valid  only  for  near- 
optimum collector  tilt  angles  and  annual  solar  radiation 
data. 

For  solar  heating  applications,  where  monthly  radia- 
tion figures  must  be  used,  the  correction  factor  may 
not  be  expressed  in  closed  form.  In  this  case.  K must 
be  determined  from  Table  1,  which  gives  monthly 
values  of  this  correction  factor  for  four  different  lati- 
tudes. Table  1 is  valid  only  for  optimum  collector  tilt 
angles. 

Recalling  the  example  on  page  I I,  the  annual  radia- 
tion for  Topeka,  KS.  was  found  as  follows:  the  annual 
radiation  level  on  a horizontal  surface  obtained  from 
Figure  B1  of  Appendix  B is  380  langleys/day.  The  lati- 
tude, i>L  (also  from  Figure  Bl ),  is  39.  Using  F.qs  4 and 
5,Q  can  be  found. 
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Q(,  = (>.0  X 10s  Btu/sq  ft/yr  (0.8  X I0H  kJ/m2/yr) 

As  .m  example  of  applying  die  monthly  correction 
(actor  used  for  heating-only  systems,  find  the  solar 
radiation  for  January  in  Topeka,  KS. 

From  Table  1.  K is  estimated  by  interpolation  between 
35  degrees  and  40  degrees  latitude  for  January.  Thus, 

k = :.06 

From  Figure  B2  of  Appendix  B,  the  mean  daily  hori- 
zontal radiation  in  January  is  190  langleys.  hence, 

[I  / 1(>0  langleys \ / 3 1 daysc  / 3,(s9  Blu  \ 

v ' day  / ’ mo  / ' sq  It-langley  ' 

llu=  2.2  X ID4  Btu/sq  ft/1110  (2.5  X I0(i  kJ/m2/mo) 
Thus, 

Qc  = KHU 

= 4.5  X 104  Btu/sq  ft/1110  (5.1  X I05  kJ/m2/mo). 


Table  1 

Monthly  Correction  Factors 
for  Heating-Only  Systems 
(for  collector  tilt  equal  to  latitude  plus  10  degrees) 


Month 

jt>° 

jy 

Latitude 

40 0 

45" 

Jan 

1.64 

1.84 

2.12 

2.49 

1 eb 

1.42 

1.55 

1.71 

1.93 

Mar 

1. 18 

1.25 

1.33 

1.44 

Apr 

.97 

1.00 

1.03 

1.07 

May 

.84 

.84 

.85 

.87 

J un 

.78 

.78 

.78 

.78 

Jut 

.80 

.80 

.80 

.81 

Aug 

.90 

.91 

.93 

.96 

Sep 

1.07 

1.1 1 

1.17 

1.24 

Oct 

1.29 

1.40 

1.5 1 

1.66 

Nov 

1.54 

1.70 

1.93 

2.23 

Dec 

1.70 

1.93 

2.24 

2.64 

ESTIMATING  SOLAR  SYSTEM 
PERFORMANCE 


General 

A critical  step  in  the  economic  analysis  of  .1  solai 
energy  system  is  a determination  ol  the  solai  system 
performance  to  find  out  how  much  of  the  energy  ic 
quireinent  is  met  by  solar  energy  and  how  much  is  met 
by  auxiliary  heat  for  a given  size  collectoi  ariay.  I11 
order  to  make  this  determination,  the  universal  curses 
arc  used. 

Analyzing  a given  solar  system  performance  requires 
the  following  steps: 

1.  Determine  the  annual  or  monthly  energy  require- 
ment. 0,  . as  described  on  pages  13-15. 

2.  Determine  the  annual  or  monthly  solar  radiation 
on  the  optimally  tilled  collector  array.  Q . for  the 
location  in  question,  using  the  method  described  be- 
ginning on  page  I 5. 

3.  Determine  the  solar  system  performance  rising 
the  universal  curve.  Two  methods  can  be  used  to 
accomplish  this: 

a.  II  a collector  aiea.  A . is  assumed,  then  the 

C 

annual  or  monthly  solar  system  performance  pata- 
meter.  P , can  be  found  from  tq  1 : 


OA 


[tq  I 


From  this  value  and  the  equation  for  the  appropriate 
universal  curve,  the  annual  or  monthly  fraction  of 
energy  met  by  the  solar  system,  p,  can  be  determined. 
Using  this  value  in  tqs  6 and  7 yields  the  energy  sup- 
plied by  auxiliary  heat,  QLA.  and  the  energy  supplied 
by  the  solar  system.  Qt  g. 

0LS  = P0L  [Lq  6 1 


If  the  analysis  is  for  heating-only  systems,  the 
monthly  0,  s and  QLA  are  obtained  from  tqs  6 and  7 
and  must  be  summed  to  get  the  annual  values  required 
for  determining  optimal  collector  area  and  making  the 
energy  and  economic  cost/benefit  analysis  as  described 
in  Chapter  4. 
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b.  The  universal  curve  can  also  indicate  what  collec- 
tor area  is  required  to  meet  a given  fraction,  p,  of  the 
annual  or  monthly  energy  requirement.  This  informa- 
tion is  useful  in  establishing  a starting  point  for  a solar 
system  analysis.  The  collector  area  can  be  determined 
bv  using  p to  find  Ps  from  the  universal  curve  equa- 
tions. hoi  P . the  collector  area  can  be  found  from: 

' (->c 

Again,  annual  or  monthly  s and  0LA  is  given  by 
hqs b and  7. 

The  next  sections  present  the  universal  curves  for 
the  various  systems  and  examples  of  how  to  use  them. 

Solar  Domestic  Hot  Water  Systems 

The  development  of  the  universal  curve  for  domestic 
hot  water  (Figure  2)  was  based  on  a two-tank  system 
in  which  the  solar-heated  tank  serves  as  a preheater  for 
a conventional  hot  water  generator.  According  to  the 
control  scheme  in  the  simulation,  the  solar  tank  is 
always  allowed  to  achieve  the  highest  temperature 
possible  under  the  given  conditions  of  solar  radiation. 
Auxiliary  energy  is  provided  to  the  conventional  hot 
water  tank  as  required  to  maintain  a preset  tempera- 
ture. (In  accordance  with  KTL  110-3-266,  the  actual 
measured  temperature  of  hot  water  delivered  to  the 
user  should  not  exceed  110°F  [43°C]).  A uniform 
annual  hot  water  demand  was  assumed  when  generating 
the  universal  curve  for  domestic  hot  water;  therefore, 
the  curve  should  not  be  used  if  the  system  being  con- 
sidered has  a highly  seasonal  demand.  For  seasonal 
demands,  heating-only  monthly  analysis  must  be  made. 

The  following  example  illustrates  application  of  the 
procedure  outlined  above  as  applied  to  solar  domestic 
water  heating. 

Example  I 

Assume  that  a barracks  at  Fort  Hood.  TX,  consumes 
30  gal  ( I 1 4 C)  of  hot  water  per  day  pel  person  and  that 
the  average  occupancy  for  the  year  is  200  people.  The 
collector  area  which  will  meet  SO  percent  of  the  load 
assuming  that  the  city  water  supply  is  55°F  (13°C)  and 
that  the  water  must  be  heated  to  1 10°F  (43°C),  can  be 
determined  in  the  following  manner: 

I.  First  determine  QL.  the  annual  energy  required 
to  heat  the  water.  Eq  3 gives: 


Q.  = DwC  (T  -T  I 

L p out  in 


(_30^)  rsons,(^^) 

L \ person-day ' ' yr  ' 

Q,  = I.OX  10s*  Btu/yr (l.l  X |0H  kJ/yi) 


2.  Estimate  (J  . the  annual  solai  radiation  on  ihc 

C 

array.  From  Fqs  4 and  5. 

cos(i>  - 7 - v ) 
°f=-^L-7)C-nu 

From  Figure  BI,  the  mean  daily  horizontal  solar  radia- 
tion at  this  location  is  445  langleys/day.  Thus. 

„ _ <445  langleyjt  <3o5daysW  3.64  Btu  \ 
v ' day  / \ yr  ' ' sq  ft-langley ' 

= 6.0  X 105  Btu/sq  ft/y r (6.8  X 106  kJ/m2/yr) 

Also  from  Figure  Bl . uL  is  33  degrees  (latitude  of  Fort 
Hood,  TX).  Since  this  is  a domestic  hot  water  system, 
the  optimum  tilt  angle.  i>c,  for  hot  water  is  equal  to 
(33  degrees).  Thus, 

Qc  =6.6  X I05  Btu/sq  ft/yr  (7.5  X 106  kJ/m2/yr). 

3.  Use  the  universal  curve  to  estimate  the  solar  per- 
formance. To  meet  50  percent  of  the  load,  p is  .5. 
From  Eq  2 (for  p equal  to  .5)  Ps  is  .88.  The  required 
collector  area  can  now  be  found  from  Eq  I . 

_ PsQl  _ (.88X1.0  X 109  Btu/yr) 
c Qc  (6.6  X 105  Btu/sq  ft-yr) 

= 1 .3  X 103  sq  ft  (120  m2) 

The  area  required  to  meet  50  percent  of  the  annual 
domestic  hot  water  ioad  is  1300  sq  ft  ( 120  m2). 

Solar  Heating  and  Cooling  Systems 

Figure  3 shows  the  universal  curve  for  solar  heating 
and  cooling  systems.  Here  the  relationship  between  Ps 
and  p is  given  by 

p = 1 - e~  32P1  07  [Eq  8J 

Analysis  of  the  system's  performance  proceeds  exactly 
as  described  for  solar  domestic  hot  water.  The  follow- 
ing example  illustrates  an  application  of  the  method 
for  solar  heating  and  cooling. 


d U.9H3N3  dVIOS  A8 
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Example  2 


The  following  example  is  for  an  administrative  build- 
ing housing  300  people  in  Champaign,  IL  (latitude  40 
degrees).  The  energy  consumption  data  were  taken 
from  monthly  utility  hillings.  The  steps  in  the  analysis 
are  ;ls  follows: 

I l ust  determine  0,  . the  annual  energy  icqiiue- 
menl  loi  heating  and  cooling. 

a.  l ot  healing,  power  records  indicate  that  for  the 
months  of  October  through  June,  144  X I08  Btu 
(152  X I08  kJ)  of  gas  were  consumed  for  heating  and 
domestic  hot  water.  When  multiplied  by  a boiler  effi- 
ciency of  .75,  the  building  heating  load  is  found  to  be 
108  X I08  Btu  (1 14  X IO8  kJ). 

b.  For  cooling,  utility  records  show  that  for  the  For  Ps  = 1 .02.  Eq  8 gives  p = .28. 
months  of  April  through  November.  82  X 104  kWh 

(28  X 108  Btu  [29.5  X 108  kJ|)  were  consumed  by  Roughly  28  percent  of  the  building's  annual  energy 

the  centrifugal  chiller.  The  chiller  load  was  found  by  requirement  for  heating  and  cooling  can  be  supplied  by 

subtracting  the  average  winter  electrical  load  from  the  the  sun.  using  50.000  sq  ft  (4645  m2)  of  collectors, 
average  summer  electrical  load  to  account  for  other, 

noncooling  loads.  Thus,  chiller  electrical  load,  when  Solar  Heating  Systems 

multiplied  by  the  COP  of  the  centrifugal  chiller  ( found  The  universal  curve  for  solar  heating  systems  ( Fig- 

froni  manufacturer’s  literature  to  be  4).  gives  a cooling  ure  4)  is  given  by  the  expression 

load  for  the  building  of  1 12  X 1 08  Btu  ( 1 18  X l08kJ).  92 

Dividing  by  the  COP  of  the  absorption  chiller  (.65)  P=  I -e~46Pi  [Eq9| 

indicates  that  1 72  X 1 08  Btu  (181  X IO8  kJ)  must  be 

supplied  to  the  absorption  chiller  to  meet  the  annual  Tq  9 is  applied  by  using  a method  similar  to  those  for 

cooling  energy  demand.  the  previous  two  systems,  except  tiiat  the  analysis  is 

carried  out  month  by  month.  Here,  the  monthly  re- 
Hence,  suits  are  summed  to  give  a seasonal  solar  percent. 

Since  the  domestic  hot  water  load  for  a heating  sys- 
tem is  small  in  comparison  to  the  space  heating  require- 
ments. it  can  be  assumed  that  the  total  hot  water  load 
is  met  by  solar  energy  during  the  nonheating  season. 
This  somewhat  simplifies  the  calculations,  in  that  the 
monthly  analysis  needs  to  be  performed  only  for  those 
months  with  significant  space  heating  loads.  The  hot 
water  energy  requirements  for  the  other  months  can  be 
summed  and  added  directly  to  QLS.  The  following 
example  shows  how  this  monthly  analysis  is  accom- 
plished for  heating-only  systems. 

Example  3 

The  administrative  building  of  Example  2 was  used 
in  this  example  to  determine  what  percent  of  the  heat- 
ing load  could  be  met  by  the  50,000  sq  ft  (4645  m2) 


QL  = 1 08  X 1 08  Btu/yr  + 1 72  X 1 08  Btu/yr 

Ql  = 280  X IO8  Btu/yr  (295  X IO8  kJ/yr). 

2.  Estimate  the  annual  solar  radiation,  Qc.  on  the 
optimally  tilted  collector  array.  From  Fqs  4 and  5, 

cosfu,  — 7 — u ) 

Q = — H 

c cos(ul  -7)  v 

Interpolating  from  Figure  Bl,  the  mean  daily  annual 
solar  radiation  for  Champaign,  IL,  is  355  langleys/day. 
Thus. 

U - /Ts-5  langleys\  / 365  daysi  / 3,69  Btu  \ 
v * day  / ' yr  ' ' sq  ft-langley  I 

= 4.8  X |05  Btu/sq  ft/yr (5.5  X 106  kJ/m2/yr) 


5 Dctciuuuc  the  solai  system  pcilniuiancc  I tom  the 
equation  Ini  the  univeisal  curve  lot  healing  and  cooling 
( I q 8).  assuming  a collector  area  ol  50.000  sq  ft  (4o45 
m2 ).  From  Eq  I . P.  is  given  by 

p ^ Ac°c  (5  X IO4  sq  ft)(5.7  X IO5  Btu/sq  ft/yr) 
s Q.  (280  X 10"  Btu/yr) 


P = 1 .02 


since  vL  = 40  degreesdhe  latitude  of  Champaign),  and 
t>c  = uL  - 10  degrees  = 30  degrees  (optimum  collector 
tilt  angle  for  heating  and  cooling). 

Thus.  0r  = 5.7  X IO5  Btu/sq  ft/yr 

(0.4  X IO1’  kJ/in2/vi ) 
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MONTHLY  INCIDENT  RADIATION /MONTHLY  ENERGY  REQUIREMENTS, 


of  collectors.  Again  the  energy  consumption  data  were 
taken  from  monthly  utilities  bills.  The  steps  in  the 
analysis  are  as  follows: 

1 . Determine  the  monthly  loads.  QL,  for  the  heating 
season.  Table  2 presents  the  building’s  natural  gas  con- 
sumption indicated  by  the  monthly  utility  records.  The 
values  of  QL  in  the  table  have  been  corrected  for  a 
boiler  efficiency  assumed  as  .75.  Note  that  the  108  X 
1118  litn/y r ( I 14  X 108  kJ/yr)  total  load  is  the  same  as 
found  in  Example  2. 

2.  Estimate  Qc,  the  annual  solar  radiation  on  the 
optimally  lilted  collector  array.  From  Eq  4. 

Uc  = K"u 

where  Hy  and  Qc  = monthly  radiation  values 

K = the  correction  factor  from  Table  1. 

Hy  is  obtained  from  the  figures  in  Appendix  B.  Monthly 
values  for  Q are  listed  in  Table  3. 

c 

3.  Determine  the  solar  system  performance  from  the 
universal  curve  for  heating  only  (Figure  4),  assuming 
a collector  area.  Ac,  of  50,000  sq  ft  (4645  m2).  From 
Eq  I.  P for  each  month  can  be  calculated  from  the 
monthly  QL  and  Qc  presented  in  Tables  2 and  3. 

From  these  values  of  P , the  monthly  values  of  p 
can  then  be  obtained  from  the  equation  for  the  univer- 
sal curve.  Finally,  from  Eq  6,  the  load  met  by  the  solar 
system,  QLS,  can  be  found  for  each  month  (Table  4). 


The  figures  in  Table  4 show  that  53.7  X 108  Btu 
(51.6  X 108  kJ)  of  heating  is  supplied  by  the  solar 
system  during  October  through  June.  The  total  heat- 
ing Ql  for  October  through  June  (from  part  1 of  this 
example)  is  108  X 108  Btu  (1 14  X 108  kJ).  Hence,  the 
solar  system  provides  50  percent  of  the  annual  heating 
requirement. 

For  the  3 months  not  considered  above.it  is  assumed 
that  all  domestic  hot  water  heating  is  performed  by  the 
solar  energy  system.  Thus, 

Qus  = 53.7  X I08  Btu  + 3(.(>  X I08  Btu) 

= 55.5  X I08  Btu  (63.5  X I08  kJ) 
and 

Ql  = 107.9  X 108  Btu  + 3(  6 X 108  Btu) 

= 109.7  X 108  Btu  (1 15.7  X 108  kJ). 

Hence,  the  annual  p is 


Thus,  the  50,000  sq  ft  (4645  m2)  of  collectors  will 
meet  51  percent  of  the  heating  and  hot  water  load  for 
a total  energy  savings  of  56  X 108  Btu  (53.4  X 108  kJ) 
compared  with  Example  2,  where  the  same  50.000  sq 
ft  (4645  m2)  of  collectors  met  only  28  percent  of  the 
heating  and  cooling  load. 


Table  2 

Natural  Gas  Consumption  for  Example  Administration  Building 


Month 

0L  (Heating) 

q,  (»:„iv, 

,'ster) 

Total  Ql 

Btu  X 10* 

Btu  X 10* 

Btu/mo  x 10* 

(kJ  x 10*) 

(kJ  x 10*) 

(kJ/mo) 

Oct 

2.0 

(2.1) 

.6 

(.6) 

2.6 

(2.7) 

Nov 

7.8 

(8.2) 

.6 

(.6) 

8.4 

(8.9) 

Dec 

12.3 

(13.0) 

.6 

(.6) 

12.9 

(13.6) 

Jan 

19.9 

(21.0) 

.6 

(.6) 

20.5 

(21.6) 

Feb 

18.5 

(19.5) 

.6 

(.6) 

19.1 

(20.2) 

Mar 

16.3 

(17.2) 

.6 

(.6) 

16.9 

(17.8) 

Apr 

12.9 

(13.6) 

.6 

(.6) 

13.5 

(14.2) 

May 

8.8 

(9.3) 

.6 

(.6) 

9.4 

(9.9) 

Jun 

4.0 

(4.2) 

.6 

(.6) 

4.6 

(4.9) 

102.5  X 

10*  Btu/yr 

107.9  X 10*  Btu/yr 

(108.1  X 10'  kJ/yr) 

(113.8  X 

10‘  kJ/yr) 
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Table  3 

Monthly  Values  for  Q(. 


Month 

l.angley /(I  ay 

B(u  sq  ft  /mo 
(kJ/m’/mo) 

k 

Dlusq  II  mo 

(kj/m’/mo) 

Oct 

275 

3.1  x 10* 
(3.5  x 10') 

1.51 

4 7 X 10* 

(5.3  X 10’ ) 

Nov 

175 

1.9  x III* 
(2.2  x It)' ) 

1.93 

3.7  x 10* 

(4.2  X 10' ) 

Dec 

135 

1.5  X ID* 
(1.7  '■  10') 

2.24 

3.4  x 10* 

(3.8  x 10’) 

Jail 

155 

1.8  X 10* 
(2.0  x 10') 

2.12 

3.8  X 10* 

(4.3  X 10’ ) 

1 eb 

240 

2.5  x 10* 
(2.8  x 10') 

1.71 

4.3  X 10* 

<4.9  x 10’ ) 

Mai 

330 

3.8  x 10* 
<4.3  x 10s ) 

1.33 

5.1  x 10* 

(5.7  x 10’ ) 

Apr 

400 

4.4  x 10* 
(5.0  x 10') 

1.03 

4.5  X 10* 

(5.1  X 10s) 

May 

510 

5.8  x 10* 
<6.6  x 10' ) 

.85 

4.9  x 10* 

(5.6  X 10’  1 

Jun 

550 

6.1  X 10* 
<6.9  X 10’) 

.78 

4.8  X 10* 

(5.4  X 10’ ) 

Table  4 

Monthly  Values  of  e and  Q 

LS 

Month 

<?L 

Btu/ino 

(kj/mo) 

Qc 

btu/sq  ft/mo 
(kj/m’/mo) 

ps 

e 

^LS 

Btu 

(k.l) 

Oct 

2.6  X 10" 

(2.7  X 10") 

4.7  X 10* 

(5.3  X 10’ ) 

9.0 

.97 

2.5  x 10" 
(2.6  X 10") 

Nov 

8.4  X 10' 

(8.9  X 10') 

3.7  X 10* 

(4.3  x 10') 

2.2 

.61 

5.1  X 10" 
(5.4  x 10" ) 

Dec 

12.9  X 10" 

( 1 3.6  X 10' 1 

3.4  x 10* 

(4.0  X 10’) 

1.3 

.44 

S.7  X 10* 
(6.0  X 10"  1 

Jan 

20.5  X 10' 

(21.6  x 10") 

3.8  X 10* 

(4.5  X 10') 

.93 

.35 

7.2  X 10" 
(7.6  X 10") 

1 eb 

19.1  x 10" 

(20.2  X 10") 

4.3  X 10* 

(5.1  X 10’) 

1.1 

.39 

7.4  X 10* 
(7.8  X 10") 

Mar 

16.9  X 10" 

(17.8  X 10") 

5.1  X 10* 

(5.6  X 10’) 

1.5 

.49 

8.3  X 10* 
(8.7  X 10" ) 

Apr 

13.5  X 10" 

(14.2  X 10") 

4.5  X 10* 

(5.0  X 10’) 

1.7 

.53 

7.2  X 10" 
(7.6  X 10"  1 

May 

9.4  X 10" 

(9.9  X 10") 

4.9  X 10* 

(5.6  X 10’) 

2.6 

.67 

6.3  X 10* 
(6.6  X 1(1"  1 

Jun 

4.6  X 10" 

(4.9  X 10") 

4.8  X 10* 

(5.5  x 10’) 

5.2 

.87 

4.0  X 10" 
(4.2  X 10") 

Total  0|-s 

53.7  X 10*  Btu 
(57  X 10" kJ) 

[ • I ■ 
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Corrections  for  Different  Collector  Types 

As  previously  described,  the  univei sal  curves  were 
developed  loi  the  leleieuec  solar  collector.  For  collec- 
tors with  different  absorber  plate  characteristics  and 
different  numbers  of  covers,  fables  5 and  6 can  be 
used  to  correct  the  results  of  the  universal  curve.  These 
tables  give  values  for  a collector  area  multiplication 
factor  which  can  be  applied  when  a particular  collector 
other  than  the  reference  collector  is  used.  In  these 
tables,  « is  the  absorptivity  and  e is  the  emissivity  of 
the  absorber  plate.  The  transmittance  of  each  cover 
was  assumed  to  be  0.9,  so  that  for  two  covers,  the 
transmittance  equals  0.81 . 

These  tables  are  useful  because  they  permit  applica- 
tion of  the  performance  curves  for  the  reference  flat- 
plate  collector  (a  single-cover  collector  with  an  absorp- 
tivity of  0.9  and  an  emissivity  of  0.10)  to  other  flat- 
plate  collectors.  For  example,  in  order  to  achieve  equi- 
valent performance  for  heating  and  cooling  from  a 
two-cover  (N  = 2),  nonselective  collector  (a  = e=  .96), 
consulting  Table  5,  one  would  multiply  the  appropriate 
collector  area  for  the  reference  collector  by  1.09  or 
increase  the  collector  area  by  9 percent. 


Table  5 

Collector  Area  Multiplying  Factor 
for  Different  Collector  Designs 
for  Heating  and  Cooling  Systems 


a 

0.96 

0.94 

0.90 

€ 

0.96 

0.30 

0.10 

1 

1.55 

1.09 

1 

2 

1.09 

0.97 

0.93 

a - Absorptivity 

e = Fmissivity 

N = Number  of  glass  covers 


Table  6 

Collector  Area  Multiplying  Factor 
for  Different  Collector  Designs 
for  Domestic  Hot  Water  Heating 
and  Heating-Only  Systems 


a 

.96 

.90 

c 

.96 

.10 

1 

1.26 

1 

2 

.96 

.93 

a - Absorptivity 

e = F.missivity 

N = Number  of  glass  covers 


A DETERMINING  THE  OPTIMAL 
^ COLLECTOR  AREA 


This  chapter  describes  a method  for  determining  the 
collector  area  with  the  maximum  economic  payback 
over  the  life  of  the  facility  for  which  it  is  being  consi- 
dered. The  method  is  based  on  the  computation  of  a 
benefit-to-cost  ratio  (B/C)  for  the  solar  energy  system, 
the  performance  of  which  is  determined  from  the 
appropriate  universal  curve.  The  benefit  associated 
with  the  incorporation  of  solar  energy  is  the  fuel  sav- 
ings; the  cost  is  the  increased  first  cost  of  the  solar 
energy  system  itself.  If  the  benefits  outweigh  the  costs, 
B/C  is  greater  than  1,  and  the  system  is  economically 
feasible. 

Basic  standards  and  guidelines  for  the  conduct  of  all 
economic  studies  by  and  for  the  Department  of  the 
Army  are  contained  in  AR  1 1-28. 11  Supplemental 
guidelines  for  performing  an  economic  analysis  of 
energy -related  systems  may  be  obtained  in  the  ECIP 
guidelines.12  Further  clarifications  of  the  basic  criteria 
may  be  obtained  through  normal  channels  from  the 
Office  of  the  Chief  of  Engineers.  Department  of  the 
Army  (ATTN:  DAEN-MPE-T). 

The  first  step  in  determining  the  optimal  collector 
area  is  to  estimate  the  initial  capital  costs  (1C)  of  the 
solar  energy  system  components  in  present-year  dollars. 
Only  the  cost  of  components  that  arc  not  normally 
part  of  the  conventional  heating  and  cooling  system 
should  be  included.  For  example,  the  cost  of  the  build- 
ing’s air-handling  system  would  not  be  considered,  but 
the  difference  between  the  cost  of  a more  expensive 
absorption  chiller  and  a lessexpensive  centrifugal  chiller 
should  be  charged  to  the  solar  system. 

Certain  cost  elements  for  solar  energy  systems  vary 
according  to  the  size  of  the  system,  while  others  are 
relatively  fixed,  regardless  of  collector  area  or  tank  size. 
The  costs  of  the  collector  and  storage  tank  are  obvious 
examples  of  items  which  are  dependent  on  collector 
area.  (The  tank  volume  is  assumed  to  be  16  lb  of  water 


11  Economic  Analysis  and  Program  Evaluation  for  Resource 
Management , AR  It -28  (Department  of  the  Army,  2 December 
1975). 

13 Letter  from  David  M.  Crabtree,  Chief,  Util.  F.ngr.  and  Ops. 
Div.,  Facilities  Engineering  (DAEN-FFU),  Subject:  Energy 
Conservation  Investment  Program  (FC1H)  Guidance  (7  Novem- 
ber 1977). 
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per  square  fool  of  collector  1 78  kg/m2 1 .)  Other  exam- 
ples include  heat  exchanger,  pump,  and  piping  costs. 
Central  system  costs  associated  with  the  solai  energy 
system,  however,  are  largely  independent  ol  collector 
area.  I he  cost  ditterenee  associated  with  the  purchase 
and  installation  of  an  absorption  chiller  is  also  relatively 
independent  of  solar  collector  area,  since  the  selection 
of  an  appropriate  absorption  chiller  is  dictated  by  peak 
building  cooling  load  for  all  but  the  smallest  collectoi 
areas.  It  is  important  that  these  costs  be  apportioned 
accurately,  since  the  benefit-to-cost  ratio  is  computed 
as  a (unction  ol  the  area  of  the  solar  collector  array. 

When  determining  the  capital  cost  of  the  solar  sys- 
tem, the  user  must  specify  the  type  of  auxiliary  energy 
system  to  be  used.  For  example,  most  solar  systems 
have  a complete  conventional  system  as  an  auxiliary. 
For  solar  systems  containing  an  absorption  chiller,  the 
auxiliary  source  of  cooling  can  either  be  a boiler  (to 
drive  the  absorption  chiller ) ora  redundant  centrifugal 
chiller.  In  the  first  case,  the  capital  cost  is  low,  since 
a credit  can  be  taken  for  the  centrifugal  chiller  and  the 
boiler  is  always  required  for  the  heating  system.  How- 
ever. the  fuel  costs  may  be  high,  since  the  absorption 
chiller  COP  is  approximately  .65.  In  the  second  case. 


cost  of  cooling  energy  (note  that  the  COP  ol  the  con- 
ventional electrical  chilling  device  must  be  considered 
when  making  this  conversion). 

,C  Divide  the  tequiied  input  heating  energy  and.oi 
domestic  hot  water  energy  by  the  heating  value  ol  the 
fuel  used  and  multiply  by  the  unit  fuel  price*  to  deter  - 
mine  the  cost  ol  heating  (note  that  in  determining  the 
heat  energy  required,  the  efficiency  of  the  boiler  01 
furnace  must  be  considered). 

4.  Convert  the  amount  of  fuel  cost  to  the  present 
worth  luel  cost  for  the  conventional  system.  PWFC  . 

C 

which  is  given  by 

PWFC  = Ml  . 

C C 

where  M is  the  cumulative,  discounted  single  amount 
factor,  and  is  expressed  by 

M = 1/2  V — ! — + ! I 

- ~0  <l-i)n  <l-i)n-r| 

Here,  n is  the  facility  life,  1 the  discount  rate,  and  t is 
the  long-term  differential  inflation  factor  for  the  fuel 
being  considered.  (Values  for  I and  t are  also  listed  m 


the  capital  costs  are  higher,  since  credit  for  the  centri- 
fugal chiller  cannot  be  taken;  however,  since  this  chiller 
operates  with  a COP  of  approximately  4,  the  fuel  costs 
may  be  considerably  lower.  Note  that  in  the  first  case, 
heating  fuel  is  the  source  of  auxiliary  energy,  while  in 
the  second  case,  the  source  is  electrical  power. 

It  the  project  under  consideration  is  programmed 
for  funding  in  the  future,  the  initial  capital  cost  (1C)  of 
the  system,  expressed  in  today's  dollars,  must  be  esca- 
lated to  the  value  of  dollars  at  that  date.  The  escalated 
1C  is  called  the  current  wording  estimate  (CWH);  the 
factors  for  converting  1C  to  CWE  are  given  in  the  ECIP 
guidance. 

Once  the  CWE  has  been  determined,  the  second  step 
in  the  economic  analysis  is  to  establish  tiie  annual 
fuel  cost  for  the  conventional  energy  system  (FC)  as 
fallows: 

1.  Determine  the  annual  heating  and/or  cooling 
energy  requirement  for  conventional  equipment  as  des- 
cribed in  Chapter  2. 

2.  Convert  the  cooling  energy  load  to  kilowatt  hours 
and  multiply  by  the  local  electrical  rale*  to  obtain  the 

•Energy  costs  here  are  fuel  prices  escalated  to  the  end  of 
the  year  in  which  the  project  is  programmed  for  funding. 


the  ECIP  guidance.)  If  more  than  one  type  of  fuel  is 
being  used,  the  PWF('c  must  be  calculated  for  each 
type,  and  the  results  summed  to  give  the  total  PWFC  . 
The  ECIP  guidance  provides  tables  of  values  of  M for 
various  fuel  types. 

The  third  step  in  the  economic  analysis  is  to  calculate 
the  present  worth  luel  cost  for  energy  required  by  the 
conventional  system  if  lire  solar  energy  equipment  is 
present,  PWFCs.  To  do  this  it  is  necessary  to: 

1.  Determine  the  total  annual  energy  requirement. 
QLA ■ for  a given  collector  area  and  systems  type  using 
the  methods  described  in  Chapter  3. 

2.  Determine  the  annual  cost  for  this  auxiliary  ener- 
gy. E\*  For  heating  and  cooling  systems,  assume  that 
solar  energy  is  used  first  to  meet  the  domestic  hot 
water  demand,  then  the  heating  load,  and  then  the 
cooling  load.  Note  in  particular  that  if  centrifugal  or 
other  electrical  chilling  devices  are  used  to  augment  the 
solar  cooling  equipment,  the  appropriate  COI's  for 
absorption  and  electrical  refrigeration  equipment  must 
be  used  to  determine  the  auxiliary  energy  that  will  be 


♦Energy  costs  are  based  on  fuel  prices  escalated  to  the  end 
ol  the  year  in  which  the  project  is  programmed  lot  funding. 
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required  lu  meet  the  cooling  demand  not  met  by  the 
solai  system.  II  a boiler  is  used  to  drive  the  absorption 
chiller,  the  boiler  efficiency  must  be  considered. 

3.  Convert  the  annual  auxiliary  fuel  cost  to  present 
worth  fuel  cost  as  before: 

PWFC  = MF 

s s 

where  M is  the  cumulative,  discounted  single  amount 
factor. 

The  final  step  in  the  economics  analysis  is  calculation 
of  the  benefit-to-cost  ratio,  B/C,  for  various  collector 
areas.  B/C  is  given  by 


PWFC  -PWFC 


where  PWFC  . PWFC  . and  CWF  are  as  defined  pre- 
viously.  Systems  are  cost-effective  only  if  this  ratio  is 
greater  than  one.  Figure  5 is  a hypothetical  plot  of  B/C 
as  a function  of  collector  area. 

The  shape  of  this  curve  is  readily  explained.  For 
very  small  collector  areas,  very  little  fuel  is  saved  to 
offset  the  fixed  costs  of  the  system,  such  as  for  con- 
trols, added  piping,  and  pumps.  As  the  collector  area 
increases,  however,  the  fuel  savings  begin  to  offset  the 
fixed  costs,  and  B/C  increases.  As  more  collectors  are 
added,  a point  of  diminishing  returns  is  reached;  addi- 
tional collector  area  provides  very  little  added  fuel 
savings.  Thus,  the  added  cost  of  these  collectors  and 
tank  cannot  be  offset  by  their  fuel  savings,  and  there- 
fore the  benefits  are  outweighed  by  the  added  costs  of 
these  components.  Hence.  B/C  begins  to  decrease. 

When  B/C  is  plotted  for  the  solar  system  being  con- 
sidered (as  in  Figure  5),  the  optimum  collector  area 
becomes  obvious— the  area  at  which  B/C  is  maximized. 
Note  that  Figure  5 concisely  summarizes  the  economics 
of  solar  systems.  Only  when  B/C  is  greater  than  1 can 
the  solar  energy  system  be  economically  justified.  In 
fact,  if  the  collector  and  tank  costs  are  too  high,  the 
fuel  savings  will  never  offset  the  fixed  costs  and  the 
B/C  will  never  achieve  a value  of  1 . 

An  application  of  the  methodology  presented  in  the 
preceding  three  chapters  enables  an  energy  and  eco- 
nomic analysis  of  the  solar  energy  system  being  consi- 
dered. If,  based  on  this  analysis,  the  solar  system  is 
found  to  be  feasible,  then  a detailed  design  analysis  can 
be  made  using  the  computer-aided  method  described  in 
Chapter  5. 


5 PERFORMING  COMPUTER  AIDED 
DESIGN  CALCULATIONS 


Introduction 

The  previous  chapter  described  a method  lor  graph- 
ically estimating  the  performance  ol  candidate  solar 
domestic  hot  water,  building  heating,  or  building  heat- 
ing and  cooling  systems.  This  method  is  sufficiently 
accurate  fora  preliminary  energy  and  economic  analysis 
of  solar  energy  systems.  For  most  projects,  however, 
selecting  the  optimum  collector  area  and  type  and 
insuring  the  selection  of  the  appropriate  collector  tilt 
angle  and  storage  tank  volume  may  require  a more 
detailed  design  analysis.  The  methods  described  in  this 
chapter  are  based  on  the  application  of  the  BLAST 
program  developed  at  CF.R1..  This  program  permits 
hourly  simulation  of  building  toads,  air-handling  sys- 
tems meeting  these  loads  and  simulation  of  central 
energy  plan  equipment,  including  solar  energy  systems, 
supplying  the  air-handlers.  Designers  who  wish  to  use 
this  program  should  obtain  the  users  and  reference 
manuals  which  describe  the  preparation  of  the  input 
required  for  using  the  program  and  the  algorithms  em- 
ployed in  performing  the  simulation. 

The  optimum  design  using  the  simulation  model  is 
again  determined  by  performing  the  economic  analysis 
as  described  in  Chapter  4.  The  simulation  tools,  how- 
ever, allow  better  determination  of  system  performance 
and,  thus,  a more  accurate  measure  of  p.  the  fraction 
of  the  load  met  by  the  solar  system  for  a given  collector 
array.  Also,  other  variables  such  as  collector  tilt  angle, 
storage  tank  volume,  and  type  of  collector  can  be  varied 
to  maximize  solar  system  performance.  In  addition,  the 
BLAST  program  allows  for  the  analysis  to  be  performed 
using  more  accurate  values  of  the  building  energy  re- 
quirements, Ql.  since  assumptions  concerning  average 
efficiency  and  COP  do  not  have  to  be  made. 

Information  Required  to  Use 
the  BLAST  Program 

To  use  BLAST,  hourly  weather  data,  including  solar 
radiation,  must  be  acquired  in  the  necessary  format  as 
input  to  the  computer  program.  While  several  versions 
(National  Weather  Service  Series  280  Solar  Data  Tapes, 
National  Weather  Service  Series  1440  Climatological 
Data  Tapes,  NBS-NOAA-ASHRAL  Test  Reference  Year 
Tapes,  and  SOLMET  Tapes)  of  suitable  weather  tapes 
exist;  obtaining  this  information  may  take  up  to  6 
weeks.  Therefore,  users  should  anticipate  use  of  the 
BLAST  program  and  request  the  tapes  early  enough  so 
they  will  be  available  when  required. 


1 lie  BLAST  program  requires  a description  of  the 
building  as  input.  The  user’s  manual  describes  the  exact 
form  of  this  input.  The  program  uses  such  information 
as  the  type  of  wall  construction,  window  area  and  con- 
struction. building  orientation,  schedule  for  occupancy, 
and  electrical  demands  in  determining  the  hourly  build- 
ing /one  loads  and  hourly  total  load  of  the  building.  If 
measured  data  are  available,  they  can  frequently  be 
used  to  check  the  building  loads  to  insure  that  para- 
meters were  properly  defined  and  reflect  the  actual 
building,  since  some  parameters,  such  as  infiltration 
and  occupancy,  aie  hard  to  determine. 

The  air-handling  systems  used  in  the  building  must 
also  be  described.  In  this  case,  such  information  as  the 
zone  air  flow  rate,  type  of  coil,  type  of  air-handling 
system,  and  type  of  control  system  used  is  specified  by 
the  user.  The  output  from  the  air  distribution  system 
simulations  is  the  hourly  energy  demand  on  the  boilers 
and  chillers  serving  the  air  handlers  in  the  building, 
v These  hourly  demands  form  the  input  to  the  central 


plant  simulation  program  which  determines  the  ex- 
pected performance  of  conventional  and  solar  energy 
system  components.  In  performing  the  central  energy 
plant  simulation,  the  user  can  specify  the  perfomtance 
of  the  solar  collector,  the  collector  area,  and  the  storage 
tank  volume  being  used  in  the  system.  The  type  of  per- 
formance data  required  for  solar  collectors  is  identical 
to  that  described  in  National  Bureau  of  Standards’ 
publication  NBSIR-74-635.13  The  slope  and  intercept 
of  solar  collector  efficiency  curves  are  the  input  para- 
meters required. 

Note  that  in  the  analysis  of  the  solar  energy  systems, 
the  BLAST  program  should  first  be  used  to  minimize 
the  energy  required  by  each  zone  by  studying  such 
effects  as  added  insulation  and  changes  in  window  area 
and  orientation.  This  should  be  followed  by  a compari- 
son of  various  air-handling  systems  and  control  schemes 
to  minimize  the  energy  demanded  by  the  boilers  and 
chillers,  assuming  these  design  changes  can  be  made  to 
the  building.  Once  this  analysis  has  been  performed, 
various  collector  systems  must  be  simulated  to  optimize 
the  performance  of  the  candidate  systems  and  to  estab- 
lish the  optimum  system  configuration. 

Once  the  load  profiles  for  the  boilers  and  chillers 
have  been  determined,  the  optimization  of  the  solar 
energy  system  can  be  initiated. 

IJJ.  S.  Hill  and  T.  Kusuda,  Methods  of  Testing  for  Rating 
Solar  Collectors  Based  on  Thermal  Performance , NBSIR-74- 
635  (National  Bureau  of  Standards,  1974). 


Setting  Up  the  Solar  Energy 
Detailed  Design  Study 

The  computer-aided  design  should  not  be  initiated 
until  the  graphical  economic  analysis  described  in  the 
previous  chapter  has  been  completed,  since  the  re- 
sults of  tl-.at  study  provide  the  starling  point  for  the 
computer-aided  study.  If  the  economic  study  has  been 
completed,  but  the  loads  determined  from  the  BLAST 
program  are  considerably  different  than  those  used  for 
analysis,  the  analysis  should  be  repeated.  The  graphical 
method  serves  to  bracket  the  area  of  consideration. 

In  order  to  optimize  the  design,  ccitain  vaiiahles  not 
previously  considered  should  be  varied  to  establish  their 
optimum  value,  as  discussed  in  the  following  section. 
The  starting  point  should  be  the  solar  system  obtained 
from  the  universal  curve  analysis  which  set  the  opti- 
mum collector  area.  The  collector  array  should  be  set 
with  an  azimuth  angle  of  180  degrees  (i.c.,  facing  due 
south)  and  tilted  at  the  optimum  angle  as  described  in 
Chapter  2.  The  tank  volume  should  be  set  at  16  Ib/sq 
ft  of  collector  (78  kg/m2). 

Methodology  for  Optimizing 
the  Solar  System  Design 

The  optimization  proceeds  in  the  following  steps: 

1.  First,  several  simulations  may  be  run  using  the 
BLAST  solar  simulation  package  with  different  tilt  and 
azimuth  angles  for  the  collector  array.  Figure  A1  illus- 
trates how  the  output  front  such  a parametric  study 
can  be  plotted  to  establish  the  optimum  tilt  for  a given 
azimuth.  A similar  plot  of  system  performance  versus 
azimuth  angle  will  indicate  the  optimum  azimuth.  This 
step  would  be  required  only  if  dictated  by  architectural 
considerations,  such  as  building  siting. 

2.  The  next  step  is  to  vary  the  collector  area  while 
keeping  the  eollector-area-to-tank-volume  ratio  constant 
(start  with  16  Ib/sq  ft  of  collector  1 78  kg/m2|)  and 
repeat  the  economic  analysis  described  in  Chapter  4. 
obtaining  a curve  similar  to  Figure  5.  The  collector  area 
with  the  greatest  B/C  ratios  is  the  optimum  collector 
area  for  this  eolleetor-area-to-tank-volume  ratio.  (Note 
that  this  is  a repeat  of  the  universal  curve  analysis  using 
the  more  precise  computer  analysis  tools.) 

3.  The  third  step  is  to  repeat  the  previous  step  lor 
different  collector-area-to-storage-tank-volumc  ratios. 
All  B/C  ratios  can  be  plotted  on  the  same  graph.  Suffi- 
cient runs  should  be  performed  to  determine  the  curve 
with  the  greatest  B/C.  If  no  other  information  is  avail- 
able, 8 and  32  Ib/sq  ft  (39  and  156  kg/m2)  are  good 
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second  ami  third  tries.  Subsequent  tries  are  based  on  the 
results  of  these  simulations.  The  curve  with  the  highest 
ratio  gives  the  optimum  tank-volume-to-collector-area 
ratio,  while  the  point  at  which  the  maximum  occurs 
establishes  the  optimum  collector  atea.  Ill  spite  of  the 
lact  tiiat  higher  values  ol  tank  insulation  will  lead  to 
largci  optimum  tank  volumes,  high  tank  insulation 
(especially  for  buried  tanks)  is  dtlltcult  to  achieve  in 
practice.  Theiefoie,  tank-voluine-to-collector  area  ratios 
which  greatly  exceed  16  Ib/sq  It  (78  xg/m‘‘)  should  be 
examined  critically  to  insure  that  the  specified  tank 
insulation  values  can  be  obtained. 

4.  The  final  step  is  to  vary  the  type  of  collector  to 
see  what  effect  this  has  on  system  life-cycle  costs.  The 
optimum  tank-volunie-to-collector-arca  ratio  found  in 
Step  d should  not  be  significantly  affected  by  the  type 
of  collector,  fire  curves  of  B/C  ratio  versus  collector 
area  should  be  constructed  as  before  by  simulating 
several  different  areas  for  each  collector  type  being 
considered. 

Note  that  the  only  differences  between  the  graphical 
and  computer-aided  methods  for  determining  system 
performance  and  optimal  collector  configuration  are 
that  the  simulation  model  provides  a more  accurate 
estimate  of  the  performance  of  a candidate  solar  energy 
system  and  permits  a more  refined  study  of  the  second- 
ary effects  of  collector  tilt,  azimuth  angle,  and  storage 
tank  capacity. 


6 PRACTICAL  CONSIDERATIONS 

Discussions  in  the  previous  chapters  focused  on  esti- 
mating performance  of  solar  collector  systems.  This 
chapter  deals  with  some  practical  considerations  in- 
volved in  ihe  final  design,  installation,  and  operation  of 
solar  energy  systems  and  is  based  on  experience  gained 
in  the  construction  and  operation  of  the  CTRL  solar 
facility. 

Freeze  Protection 

Because  the  components  in  the  collector  loop  are 
sometimes  exposed  to  conditions  of  low  outdoor  tem- 
perature. provision  must  be  made  in  many  locations  for 
freeze  protection  of  the  heat  transfer  fluid.  One  ap- 
proach to  this  problem  is  to  drain  the  collectors  each 
time  the  ambient  temperature  approaches  freezing; 
however,  in  some  climates,  particularly  for  the  larger 
systems,  this  method  may  not  be  practical.  More  often. 


a lluid  other  than  pure  water,  such  as  ethylene  glycol, 
propylene  glycol,  or  silicones,  is  circulated  in  any  ex- 
posed component.  A heat  transfer  penalty  may  he  paid, 
however,  since  these  liquids  may  have  a higher  viscosity 
and  a lower  specific  heat  than  pure  water  In  addition, 
for  domestic  hot  water  systems,  local  building  codes 
should  he  consulted  before  the  heat  traiislet  lluid  is 
chosen,  since  it  is  possible  to  contaminate  the  potable 
water  supply.  In  many  areas,  two  walls  must  separate 
any  toxic  lluid  from  potable  water  supplies 

Stagnation 

A malfunction  of  the  collector  pump  or  system  con- 
trols can  lead  to  stagnation  of  the  fluid  in  the  collection 
loop.  Under  conditions  of  high  solar  insolation,  temper- 
atures as  high  as  350°  ( 177°)  can  he  attained  in  some 
fiat  plate  collectors.  The  possibility  ol  stagnation  must 
be  anticipated  in  the  design  of  any  system,  f irst,  the 
collector  must  be  able  to  withstand  these  high  tempera- 
tures. Second,  the  collector  tin  id  should  not  decompose 
if  stagnation  occurs.  In  addition,  provision  (e.g..  an 
expansion  tank)  must  he  made  for  thermal  expansion 
of  the  collector  lluid  under  these  conditions. 

Thermal  Expansion 

Temperature  variations  in  the  collector  loop  are  often 
greater  than  those  encountered  in  normal  hydronic 
heating  systems.  This  is  because  the  collector  piping  is 
subjected  to  conditions  ranging  from  the  coldest  am- 
bient to  the  highest  fluid  temperatures.  Tor  this  reason, 
additional  caution  must  be  taken  to  insure  that  some 
type  of  mechanical  strain  relief  is  provided  to  allow  for 
differential  thermal  expansion  between  the  collector 
piping  and  ihe  collector  supports.  This  relief  com- 
monly consists  of  expansion  loops  at  intercollector 
connections. 

Flow  Balancing  in  the  Collectors 

In  some  of  the  larger  solar  systems  it  is  not  uncom- 
mon to  have  20  collectors  in  a single  array.  In  such  a 
case,  there  may  be  a significant  flow  imbalance  from 
one  collector  to  another.  Any  collectoi*  operating  at  a 
lower  than  design  flow  will  exhibit  a higher  than  design 
temperature  difference  across  them.  This,  in  turn, 
lowers  the  collection  efficiency.  Two  common  methods 
for  addressing  this  problem  are  the  incorporation  of 
flow-balancing  valves  or  the  use  of  a two-pipe  reverse- 
return  strategy. 

Venting  the  Collector  Loop 

Particular  attention  must  he  given  to  the  possibility 
of  air  accumulating  in  the  collector  array.  A frequent 
problem  encountered  in  the  past  has  been  accumula- 
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lion  of  pockets  ol  air  at  the  lop  of  the  collector  bank. 
These  air  pockets  prevent  the  fluid  from  passing  over 
the  collector  absorber  plate,  thus  reducing  the  area 
over  which  solar  energy  is  collected.  Such  “vapor  lock- 
ing" can  be  prevented  by  generous  use  ol  automatic  or 
manual  vents  in  the  collectoi  system  and  by  avoiding 
shaip  bends  01  othci  piping  arrangements  which  can 
Map  ait.  Lise  ol  an  sepaialors  in  the  collectoi  loop  is 
also  recommended . 

Prevention  of  Natural  Convection 

l ice  convection  should  be  prevented  from  dissipating 
solar  heat  through  the  collector  bank  when  the  solar 
energy  system  is  not  in  operation.  This  can  be  prevented 
by  using  two-position  valves  which  close  when  the  col- 
lector system  is  de-energized,  but  are  opened  if  the 
system  is  in  operation. 

Piping  of  Hot  Fluids 

Pumps  in  the  system  should  be  positioned  carefully, 
particularly  when  the  storage  tank  is  to  be  maintained 
at  or  near  atmospheric  pressure.  A potential  cause  of 
system  failure  is  the  cavitation  of  pumps  used  to  pump 
hot  fluids.  Piping  considerations  should  insure  that  at 
the  collection  and  storage  loops,  pumps  are  positioned 
so  that  there  is  always  a net  positive  suction.  This  will 
prevent  flashing  when  the  pumps  are  energized.  When- 
ever possible,  circulating  pumps  should  be  used  for 
closed  loops,  and  self-priming  pumps  for  open  loops. 

Storage  Tank  Insulation 

Practical  experience  indicates  that  heat  loss  from 
insulated  tanks  is  frequently  higher  than  predicted.  For 
buried  tanks  in  particular,  special  care  must  be  taken  to 
insure  that  an  impenetrable  vapor  barrier  surrounds  the 
insulation  of  the  tank.  In  addition,  storage  containers 
with  a high  volume  to  area  ratio  should  be  considered 
to  minimize  the  heat  loss.  During  the  heating  season, 
this  consideration  is  less  important  for  tanks  located 
inside  the  building,  since  any  heat  loss  from  the  storage 
tank  heats  the  mechanical  room.  However,  in  the  sum- 
mer such  a tank  may  contribute  to  the  building’s  cool- 
ing load,  in  addition  to  occupying  valuable  floor  space. 
These  factors  should  be  considered  when  designing  and 
placing  the  storage  tank. 

Piping  of  the  Auxiliary  Energy  Supply 

For  practical  solar  energy  systems,  some  auxiliary 
energy  supply  will  be  required.  Since  there  are  many 
possible  configurations,  individual  cases  must  be  exam- 
ined to  determine  whete  auxiliary  energy  can  be  derived. 
However,  any  auxiliary  supply  boiler  or  other  device 
should  be  piped  parallel  with  the  solar  energy  solar 
tank.  Figure  6 illustrates  examples  of  a correct  and  in- 


correct method  lot  piping  an  auxiliary  energy  supply. 
If  the  auxiliary  supply  is  piped  in  series  with  the  tank, 
under  conditions  of  low  tank  temperature  the  auxiliary 
supply  will  not  only  he  required  to  meet  the  heat  energy 
demands  of  the  building,  but  will  also  begin  to  add 
energy  to  the  lank.  This  is  particularly  undesirable  be- 
cause n will  requite  an  unnecessarily  large  auxiliary 
liealei  capable  ol  meeting  mslanl.meous  demand  while 
simultaneously  supplying  an  added  demand  by  deliver- 
ing energy  to  the  storage  lank  It  is  also  undesirable 
because  the  stored  eneigy  should  be  derived  only  from 
the  solar  energy  system;  application  ol  auxiliary  heat 
to  the  stored  energy  raises  the  tank  temperature,  which 
reduces  collector  efficiency  and  increases  tank  heat 
losses. 

Controls 

The  objective  of  a solar  energy  system,  whether  it  is 
used  for  hot  water,  heating,  or  cooling,  is  to  maximize 
the  collection  and  storage  of  energy  from  the  sun.  In 
this  regard,  the  importance  of  the  control  strategy  and 
equipment  cannot  be  overemphasized.  While  the  con- 
trol of  these  systems  typically  does  not  require  new 
technology,  special  care  must  be  exercised  in  determin- 
ing when  the  collector  pump  is  actuated,  because  it  is 
this  unit  which  allows  energy  to  be  transferred  to  the 
storage  tank. 

The  most  straightforward  method  for  initiating  How 
in  the  collection  loop  would  be  to  sense  the  temperature 
difference  between  the  collector  plate  surface  and  the 
storage  tank.  If  a positive  difference  existed  (i.e.,  the 
collector  is  hotter  than  the  tank),  solar  energy  would 
be  available  and  the  collector  system  could  be  energized. 
Unfortunately,  if  this  strategy  were  adopted,  the  sys- 
tem would  cycle  excessively  as  the  collector  tempera- 
ture oscillated  above  and  beiow  the  tank  temperature. 

This  problem  is  commonly  attacked  through  the  use 
of  hysteresis  in  the  turn-on,  turn-off  temperatures;  that 
is,  the  collector  pump  is  activated  when  the  collector 
plate  is  15°F  (8.3°C)  warmer  than  the  lank,  but  it  is 
de-energized  when  it  is  only  3°F  (|7°C)  warmer.  In 
this  strategy  it  is  important  that  the  control  system 
have  the  capability  for  sensing  temperature  differences 
accurately  over  the  entire  region  of  control.  When  the 
tank  is  at  180°F  (82.2°C)  during  the  summer,  for 
example,  collection  would  be  initiated  at  a plate  tem- 
perature of  195°F  (90.6°C).  During  the  winter,  when 
the  storage  is  100°F  (43.3°C).  pumping  should  begin 
when  the  plate  reaches  125°F  (51 ,7°C ).  Sensing  these 
temperature  differences  accurately  over  this  wide  oper- 
ating range  places  an  additional  burden  on  the  system 
control  electronics. 
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HEAT  EXCHANGER 


a.  CORRECT 


Figure  6.  Piping  tire  auxiliary  energy  supply. 


While  the  use  of  hysteresis  is  the  most  common  ap- 
proach to  the  problem  of  differential  temperature  con- 
trol, another  method -proportional  control-has  been 
employed.  Here,  the  speed  of  pumping  is  varied  in  pro- 
portion to  the  difference  in  plate-tank  temperature. 
The  advantage  in  this  case  is  that  energy  is  collected  on 
days  when  the  solar  insolation  is  high  enough  to  heat 
the  plate,  but  not  to  a temperature  1 5°  F (8.3°C)  greater 
than  the  storage.  The  disadvantage  lies  in  the  fact  that 
large,  variable-speed  pumps  are  not  common.  Propor- 
tional control  is  used  most  frequently  in  residential- 
scale  systems. 

Sizing  of  Auxiliary  Equipment 

As  stated  previously,  solar  energy  systems  are  nor- 
mally backed  by  auxiliary  equipment  capable  of  meet- 
ing a building's  heating  and  cooling  loads.  (This  require- 
ment ol  a redundant  HVAC  system  is  one  reason  that 
the  initial  cost  of  using  solar  energy  is  so  high.)  For 
heating-only  applications,  the  first  cost  and  part-load 
characteristics  of  boilers  are  such  that  the  economic 
penalty  paid  for  having  a full-sized  auxiliary  is  not  too 
great.  Because  this  is  not  necessarily  the  case  for  ab- 
sorption air  conditioning,  the  sizing  of  the  chiller  back- 
up system  becomes  more  critical. 


During  the  cooling  season,  there  will  be  times  when 
the  auxiliary  cooling  system,  in  supplementing  the 
solar  air-conditioner,  will  be  operating  at  a small  fraction 
of  its  full  load  capacity.  In  terms  of  energy  saved,  the 
low  efficiency  of  the  auxiliary  chiller  under  these  con- 
ditions refutes  the  performance  ot  the  overall  system. 
However,  if  the  back-up  chiller  has  been  sized  too  small, 
it  might  not  be  able  to  meet  the  load  during  extended 
periods  of  hot.  cloudy  weather. 

One  approach  to  the  problem  of  equipment  size  has 
been  to  have  both  a large  and  small  back-up  chiller. 
This  allows  both  units  to  be  operated  efficiently  over  a 
wide  range  of  cooling  loads.  (The  first  cost  would  be 
higher  here  too.)  At  the  very  least,  it  is  important  to 
realize  that  an  effort  should  be  made  to  insure  that 
auxiliary  chilling  equipment  is  not  greatly  oversized. 

Absorption  Air  Conditioning 

The  performance  of  systems  involving  solar  absorp- 
tion air-conditioning  has  sometimes  been  less  than 
expected.  At  least  two  factors  account  for  this.  The 
first  of  these,  excessive  machine  cycling,  is  essentially  a 
control  problem.  (Cycling  penalizes  the  equipment  per- 
formance by  not  allowing  the  machine  sufficient  time 


' 
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to  reach  its  full-rated  capacity.)  For  example,  if  the 
absorption  chiller  is  activated  when  the  tank  tempera- 
ture is  greater  than  185°F  (85°C),  but  de-energized  at 
temperatures  less  than  180°F  (82.2°C),  a five-degree 
hysteresis  or  dead  band  is  said  to  exist  in  this  control 
parameter.  If  the  hysteresis  is  increased  (i.e.,  the  on- 
off  temperature  difference  is  made  larger),  the  chiller 
will  come  on  less  frequently,  but  for  longer  periods  of 
time,  theiehy  reducing  the  ill  effects  of  cycling. 


7 CONCLUSIONS 

The  following  conclusions  can  be  drawn  concerning 
the  methods  presented  in  this  report: 

1.  The  graphical  method  of  evaluating  solar  energy 
systems  is  based  on  established  technology  and  calcula- 
tion methods,  and  provides  sufficient  accuracy  tin  use 
in  feasibility  studies. 


fhe  second  problem  associated  with  solar  cooling 
arises  when  the  absorption  chillei  is  underfired.  ( For  a 
given  condensing  watet  temperature,  there  is  a required 
minimum  generator  water  temperature,  below  which 
the  unit  is  said  to  be  underfired.)  The  penalty  under 
these  conditions  is  that  the  capacity  of  the  chiller  is 
reduced  greatly.  Underfiring  can  best  be  avoided  by 
altering  the  control  strategy  so  that  the  minimum  tem- 
perature of  water  pumped  to  (he  generator  is  raised. 
For  example,  the  tank  cut-off  temperature  for  cooling 
could  be  raised  from  180°F  to  1 85°F(82.2°C  to  85°C). 
Strictly  speaking,  raising  the  tank  cut-off  temperature 
is  not  the  best  approach,  since  there  may  be  a range  of 
condensing  water  temperatures  for  which  a tank  tem- 
perature of  180°F  (82.2°C)  does  not  cause  the  unit  to 
be  underfired.  A more  accurate  method  would  be  to 
have  a control  system  capable  of  measuring  condensing 
and  generating  water  temperature  simultaneously,  and 
subsequently  running  the  chiller  within  its  specified 
operating  region. 


2.  The  method  for  economic  analysis  includes  all 
cost  factors  which  must  be  considered  m evaluating 
solar  energy  systems. 

3.  The  computer  simulation  method  is  based  on 
established  calculation  methods  and  can  be  used  in  pre- 
paring design  instructions  and  evaluating  final  designs. 

4.  The  input  data  required  for  both  the  graphical 
and  computer  simulation  methods  can  be  readily  ob- 
tained by  District  personnel. 

8 RECOMMENDATIONS 

It  is  recommended  that  District  Engineers  use  the 
procedures  contained  in  this  report  for  use  in  designing 
solar  energy  systems;  in  addition,  it  is  recommended 
that  Facilities  Engineers  use  these  procedures  for  pre- 
paring Dl)  1341  submissions  on  solar  projects. 
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FRACTION  OF  ENERGY  REQUIREMENTS 
MET  BY  SOLAR  ENERGY 


APPENDIX  A: 

EFFECTS  OF  COLLECTOR  TILT,  HEAT  EXCHANGER 
EFFECTIVENESS,  TANK  VOLUME,  AND  TANK  INSULATION  ON 
THE  THEORETICAL  PERFORMANCE  OF  SOLAR  ENERGY  SYSTEMS 

In  Chapter  2 il  w:is  staled  that  the  single  most  important  factor  m determining  the 
performance  of  a solai  energy  system  is  I he  assumed  eolleeloi  efficiency.  In  support  ol 
•hat  statement,  Figures  At  through  A4  present  graphs  of  percent  solar,  p,  versus  P which 
explicitly  shows  the  dependence  ol'p  on  collector  array  tilt,  heat  exchanger  effectiveness, 
tank  volume,  and  lank  insulation. 
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Figure  Al.  Dependence  of  solar  fraction  on  collector  slope. 
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Figure  B6.  Mean  Jaily  solar  radiation  for  May.  Conversion  factor: 
I langley  = 3.64  Btu/sq  ft  or  41.84  kJ'ni2. 


onversion 


Figure  BIO.  Mean  daily  solar  radiation  for  September.  Conversion  factor 
1 langley  = 3.69  Btu/sq  ft  or  41.84  kJ/ni2. 


APPENDIX  C: 

EXAMPLE  OF  SOLAR  SYSTEM 
PRESENT  WORTH  COST  ANALYSIS 


A design  engineer  wants  to  determine  the  eolleetor 
area  which  has  the  greatest  economic  B/C  ratio14  for 
solar  heating  and  cooling  of  a new  administration  build- 
ing similar  to  the  one  used  for  Examples  2 and  3 in 
Chapter  3.  For  the  purpose  of  this  example,  he/she  is 
to  compare  the  costs  of  a building  with  conventional 
oil  heating  and  electric  centrifugal  cooling  to  a building 
heated  and  cooled  by  a solar  energy  system.  In  the 
solar  case,  auxiliary  heating  is  supplied  by  an  oil-fired 
boiler  and  auxiliary  cooling  is  provided  by  an  electric 
centrifugal  chiller.  For  this  example,  it  is  assumed  that 
the  analysis  is  being  performed  in  FY78  for  a project 
programmed  for  funding  in  FY80. 

Input  Data: 

From  Examples  2 and  3,  the  annual  loads  are 

Building  Heating  Load  10.8  X 103  mBtu/yr* 
Building  Cooling  Load  1 1.2  X 103  mBtu/yr 

The  FY78  cost  of  energy  is  assumed  to  be 

Price  of  Electricity  $.027/kWh 

Price  of  Oil  $.50/gal**+ 

The  energy  costs  escalated  to  the  end  of  FY8015  are 

Electricity  ($.027)(1.16)2  = S.036/kWh 
Oil  ($.50)(I.16)2  = $. 673/gal 

Heating  and  Cooling 
Step  1 

Following  the  procedures  in  Chapter  4,  the  capital 
costs  of  the  solar  system  components  must  be  deter- 


" l etter  from  David  M.  Crabtree.  Chief,  Util.  Engr.  and  Ops. 
I)iv„  facilities  Engineering  (DAEN-FEU),  Subject:  Energy 
Conservation  Investment  Program  (ECIP)  Guidance  (7  Novem- 
ber 1977). 

*SI  conversion  factor:  I Htu  = 1.055  kJ.  (To  avoid  confu- 
sion. SI  equivalents  are  not  given  in  (his  example;  instead,  an 
SI  conversion  factor  is  given  t he  first  time  a particular  unit  is 
used.) 

**S1  conversion  factor:  I gal  = 3.785  C. 

This  price  is  used  for  example  purposes  only;  it  may  be 
somewhat  higher  than  the  actual  fuel  cost. 

15  Letter  from  David  M.  Crabtree  (7  November  1977). 


mined.  Since  in  this  case  the  absoipnon  chillci  is  ie- 
dundant.  its  entire  pi  ice  is  charged  to  the  solai  system 
Table  Cl  summarizes  estimated  costs:  design  costs  have 
been  ignored  for  the  purposes  of  this  example. 

Since  the  building  of  Examples  2 and  3 contains  a 
385-ton  cooling  unit,  the  cost  of  an  absorption  chiller 
($20,000)  is  listed  in  the  table. 

The  cost  of  controls,  sensors,  valves,  and  minimal 
piping  is  nearly  independent  of  system  si/e  and  is  esti- 
mated at  roughly  $30,000. 

The  remaindei  of  the  items  m the  table  are  depen- 
dent on  collector  area  (A  ).  These  values  are  meant  as 

C 

estimates  tor  the  purposes  ot  the  example  only.  Thus, 
the  initial  capital  costs  (1C)  are 

IC=  5.0  X I04  + 14.18  A . 

c 

Escalating  this  cost  to  the  end  of  FY8016  gives  the 
current  working  estimate  (CWE) 

CWE  = (5  X 10“  + 14.1 8 A K1.07K  1.065) 

C 

CWE  = 5.7  X 104+  16.2  A 

C 

Step  2 

Once  the  cost  data  and  CWE  have  been  determined, 
the  second  step  in  the  comparison  procedure  is  the 
establishment  of  the  annual  fuel  cost  for  the  conven- 
tional energy  systems.  As  described  in  Chapter  4,  con- 
ventional system  annual  fuel  costs  are  estimated  as 
follows: 

1 . Determine  the  annual  heating  and/or  cooling 
energy  requirement  for  conventional  equipment  as  des- 
cribed in  Chapter  2. 

For  this  example  the  heating  load  is  10.8  X I03 
mBtu/yr  and  the  cooling  load  is  1 1.2  X 103  mBtu/yr 
(shown  in  Examples  2 and  3 of  Chapter  3). 

2.  Convert  the  cooling  energy  load  to  kilowatt  hours 
and  multiply  by  the  local  electrical  rate  to  obtain  the 
cost  of  cooling  energy  (note  that  the  coefficient  of  per- 
formance of  the  conventional  electrical  chilling  device 
must  be  considered  in  making  this  conversion). 


Letter  from  David  M.  Crabtree.  Chief,  Util.  Engr.  and  Ops. 
Div.,  facilities  Engineering  (DAEN-FEU),  Subject  Energy 
Conservation  Investment  Program  (ECIP)  Guidance  (7  Novem- 
ber 1977). 
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Tabled 

Estimated  FY78  Costs  for  Heating  and  Cooling  Example 


Absorption  Cooling  System 
C'ontrols,  Sensors,  Valves,  Piping 
Heat  Lxchanger 
Fluid 

Storage  lank 

Additional  Plumbing 

Pumps 

Collector 

Labor 


l ived  Costs 

$ 


5.0  X I04 


*SI  C onversion  lactor:  I sq  It  0.0929  in2. 


Cost/sq  ft* 
of  Collector,  $ 


For  a cooling  load  of  1 1.2  X 103  inBtu/yr  and  a 
chiller  COP  ol  4,  the  FY80  cost  of  cooling  energy  (at 
$. 036/kWh)  is: 


(1 1.2  X 103  mlilu/yrK  1/4) 


:.412X  10“'3  mBtu 


\ / S.036  \ 

But  ' ' kWh  ' 


= 29.5  X 1 03/yr 


3.  Divide  the  input  heating  energy  and/or  domestic 
hot  water  heating  required  by  the  heating  value  of  the 
fuel  used  and  multiply  by  the  unit  fuel  price  to  deter- 
mine the  cost  of  heating  (note  again  that  in  determining 
the  heat  energy  required,  the  efficiency  of  the  boiler  or 
furnace  must  be  considered). 

For  a boiler  which  is  75  percent  efficient,  using  a 
heating  value  of  .138  mBtu17/gal  for  oil,  the  FY80  cost 
of  heating  (at  S.673/gal  for  oil)  is: 

= 70.2  X 103/yr 

4.  Convert  the  annual  fuel  cost  to  present  worth 
fuel  cost  (PWFC).  This  value  provides  the  baseline  for 
comparison  of  conventional  and  solar  energy  systems. 


Let  ter  from  Davkl  M.  Crabtree.  Chief,  Util,  lingr.  and  Ops. 
Div..  facilities  engineering  (DAF.N-I  KU),  Subject:  Energy 
Conservation  Investment  Program  (KCIP)  Guidance  (7  Novem- 
ber 1977). 


The  PWFC  for  the  conventional  system  (PWFC  ) is 
given  by 

PWFC  = MF 

c c 

where  M is  the  single  amount  factoi  for  energy  prices 
over  the  life  of  the  facility.  If  the  facility  life  is  25 
years,  M is  19.72  for  oil18  and  17.67  for  electricity. 1H 
Hence,  for  the  cooling  system. 

PWFC  (elec)  = ( 1 7.67X29.5  X 103/yr) 


= .52  X 106 


and  for  heating. 


PWFC  (oil)  = ( 19.72X70.2  X 103/yr) 

= 1.38  X 106 

The  total  PWFC  is  the  sum  of  these  two. 

C 

PWFC  = 1.9  X 106 

c 

Step  3 

The  third  step  is  to  calculate  the  amount  of  auxiliary 
fuel  or  electrical  energy  required  annually  by  the  solar 
system  for  various  collector  array  sizes  based  on  the 
methods  described  in  Chapter  3.  Procedures  for  deter- 
mining the  present  wortli  of  auxiliary  fuel  cost  for 
various  solar  energy  systems  are  as  follows: 

'"Letter  from  David  M.  Crabtree  (7  November  1977). 
'’Letter  from  David  M.  Crabtree  (7  November  1977). 


1.  Determine  the  total  annual  energy  supplied  by 
the  solar  energy  system,  Q1S.  hy  the  method  prescribed 
in  Chapter  3.  Following  Example  2 of  Chapter  3,  is 
28  X I03  mBtu/yr  and  Qc  is  .57  mBtu/sq  ft/yr.  For  a 
collector  area.  A^.  of  5 X 104  sq  ft  (a  guess),  Eq  1 may 
be  used  to  find  F 


F = 1.02 

S 

For  ibis  F . the  equation  for  the  heating  and  cooling 
universal  curve  gives  a p of  .28.  Hence. 

U,  s = pQl  = .28(28  X 103  mBtu) 

Q[  S = 7.8  X I03  mBtu/yr 

This  is  the  amount  of  energy  supplied  by  the  solar 
energy  system. 

2.  Determine  the  annual  energy  cost.  F .for  auxiliary 
fuel. 

As  stated  in  Chapter  4.  solar  energy  is  used  first  to 
satisfy  the  building  heating  load,  then  the  cooling  load. 
In  order  to  calculate  the  annual  energy  cost  for  auxiliary 
fuel,  the  total  annual  building  thermal  demand.  QL. 
must  be  separated  into  the  thermal  requirements  for 
heating  and  cooling,  respectively,  and  compared  with 
the  amount  of  solar  energy,  QLg.  supplied. 

Following  Example  2 of  Chapter  3,  the  thermal 
energy  requirements  for  heating  and  cooling  for  this 
case  are  10.8  X 103  mBtu  and  17.2  X 103  mBtu, 
respectively.  (Note  that  these  sum  to  the  total  require- 
ment, Ql.)  For  a collector  area  of  5 X 104  sq  ft,  QLS 
has  been  calculated  to  be  7.8  X 103  mBtu.  This  energy 
is  first  used  to  meet  the  heating  requirement.  Because 
Qbs  is  less  than  the  heating  demand,  a total  of  (10.8 
- 7.8)  X !03  mBtu  must  be  purchased  for  heating 
auxiliary  and  I 7.2  X 103  mBtu  for  cooling. 

The  annual  cost  for  each  fuel  type  may  now  be 
computed.  For  heating. 


/ 3 X I03  in B t u\ 

1 $.673  \ 

/ 1 gal  \ ( 1 \ 

\ year  ' 

' gal  ) 

'.138  mBtu M. 75/ 

Where  .75  is  the  estimated  boiler  efficiency 
Fs (oil)  = $1.95  X 1 04/year 

Because  a centrifugal  chiller  (COP  of  4)  supplies  the 
cooling  backup,  the  cost  of  auxiliary  fuel  for  cooling  is 


F 

S 


, , . (17.2  X 10 

(elec)  = 1 


mBtu 


year 


/ S.03O  \ / kWh 

' kWh  ' ' 3.412  X 10~3  mBtu 

Ft(elec)=  $2.05  X l04/yr 


The  factor  of  .65  comes  from  the  fact  that  the  thermal 
demand  for  cooling  ( 1 7.2  X I03  mBtu)  was  computed 
on  the  basis  of  absorption  air  conditioning.  Now  that 
the  backup  system  is  assumed  to  be  conventional,  this 
factor  is  regained. 


3.  Convert  the  annual  fuel  cost  to  present  worth 
fuel  cost,  FWFC  . As  in  the  conventional  case, 

S 

PWECs  (oil)  = MF's  (oil) 

= (19.72X1.95  X 104) 

FWFC  = S3.85  X 105 

S 

and 


PWFC  (elec)  = MF  (elec) 

S'  S 

= ( 1 7.67)(2.95  X 104) 
= $5.21  X 105 


The  total  PWFCs.  for  a collector  area  of  5 X 1 04  sq  ft. 
is  the  sum  of: 


PWFC  = $3.85  X 105  + $5.21  X 105 

S 

PWFC  = 9.06  X I05 

S 


Step  4 

The  total  B/C  ratio  for  the  solar  heating  and  cooling 
system  is  given  by 

PWFC  -PWFC 

n ' S 


The  solar  energy  system  represents  a savings  over 
the  conventional  system  only  if  B/C  is  greater  than  1.0. 

From  Step  1,  for  an  Ac  of  5 X 104  sq  ft 

CWE  = 5.7  X 104  + 16.2(5  X 1 04 ) 

CWE  = 8.67  X I05 

Using  PWFC^  and  PWFC's  from  Steps  2 and  3.  B/C  can 
now  be  computed. 


1 


« 
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Ql  g = 14  0 X 10*  in  Bt  u 


Since  this  amount  of  energy  exceeds  the  heating  de- 
mand. all  backup  must  be  supplied  for  cooling.  Thus, 

F = (elec)  = ( 14  X 1 03 ) (^f)  1 U3<1 


4 ' V 3 .4 1 2 X KTV 

= $24.0  X 103. 
t his  implies  that 

I'WFC's  = (I7.67)(24.UX  I03)  = $4.2  X 10s 

For  an  A(,  of  IX  I05  sq  Cl.  Also, 

(AVI  = 5.7  X It)4  + I6.2(  I X 10*)=  1.68  X 10°. 

Using  the  same  PWF(’c  as  before. 

1.9  X I06  -4.2  X I05 


B/C  = 


B/C  = .88 


1.68  X 10 


Iff" 


Cot  a collector  area  of  100.000  sq  ft.  The  larger  collec- 
tor area  has  resulted  in  a lower  B/C  ratio. 

If  a collector  area  of  25.000  sq  ft  is  considered.  B/C 
is  computed  to  be  1.11.  Since  a maximum  has  been 
found  in  B/C,  a more  exact  analysis  of  the  problem  is 
indicated.  Hence,  the  computer  simulation  program 
should  probably  be  consulted. 

Heating  Only 

Fxperience  indicates  that  a high  percentage  of  the 
auxiliary  fuel  requirement  is  normally  used  for  cooling. 
This  seems  to  suggest  that  a heating-only  system  might 
(from  an  economic  standpoint)  be  more  feasible.  Thus, 
the  previous  example  is  repeated  considering  only  the 
heating  load. 


Step  1 

Determine  capital  costs.  Table  C2  summarizes  esti- 
mated costs.  For  this  case,  the  absorption  chillei  cost  is 
saved,  and  the  current  working  estimate.  (AVI  (IV  NO), 
is  given  by 

(AVI  = $3.42  X |04  r $I0.4<>  A_ 

Step  2 

Hie  computation  of  conventional  fuel  costs  is  based 
on  the  heating  load  alone.  For  a building  heating  load 
of  10.8  X I0J  niBiu/yr.  with  a .75  boiler  efficiency, 
the  FY80  annual  fuel  costs,  as  in  the  last  example,  are, 

Fc  = (10.8X  IU-1  mBtu/yrX  I /-75  KS. 673/gal) 

( 1 gal/.  1 38  mBtu) 

F = S70.2  X l03/yr 

C 

Hence. 

I’WFC,  = (19.72)(70.2  X 103) 

I'WI  C - 1.38  \ I O'1 

Step  3 

To  compute  the  cost  of  auxiliary  fuel  lot  the  solai 
system.  A must  be  determined  from  monthh  (J 
and  Q data.  The  values  I'm  Q,  and  Q for  the  heating 
months  are  given  in  Hxample  3.  Chapter  3.  This  infor- 
mation (along  with  an  assumed  collector  area)  allows 
calculation  of  a 1’  for  each  month.  1 he  universal  curve 

S 

for  heating  (Figure  4)  then  gives  the  p for  each  IH 
Finally  a monthly  0LS  can  be  figured  ( pUj  ) and 
summed  to  get  a yearly  Q( 

Table  C3  summarizes  the  procedure  using  the  Uj 
and  Qc  from  Example  3 of  Chapter  3.  with  a collector 
area  of  4 X I04.  The  QLA  can  be  determined  from 

QLA  =QL-Q,S=  |0.8X  I03  - 4.69  X 103  mBtu/yr 

Qla  = 6.2  X 103  mBtu/yr 

The  solar  auxiliary  fuel  cost.  F . is  given  (assuming  a 
boiler  efficiency  of  .75)  by 

Fs  = (6.2  X 103  mBtu/yrXl/-75)(S.673/gal) 

( I gal/.  1 38  mBtu) 

F = S40.3  X 103 


Table  C2 

Estimated  (FY781  Costs  tor  Heating-Only  Example 

ll  \ t*ll  C'osls  Cost  Sq  ft 

$ of  Collector.  % 

Controls.  Valves,  Sensors.  Piping  3 X I04 


Heat  Exchanger  .5 

Fluid  .18 

Tank  1.0 

Additional  Plumbing  .75 

Pumps  .75 

Collector  5.0 

Labor  1.0 


3 X I04  59.18  Ac 


Table  C3 

Summary  of  Procedures  Used  to  Determine 


Month 

Ol 

MBlu/mo 

Btu/sq  It/mo 

e 

MBtu/mo 

Oct 

2.6  x 102 

4.7  X 10" 

7.23 

94 

244 

Nov 

8.4  x I02 

3.7  X 10* 

1.76 

.54 

454 

t>ec 

12.9  X t0! 

3.4  X 10* 

1.05 

.38 

490 

Jan 

20.5  X 102 

3.8  X 10* 

.74 

.29 

595 

Feb 

19.1  x Id2 

4.3  X 10" 

.90 

.34 

649 

Mar 

16.9  x 102 

5.1  X 10" 

1.21 

.42 

710 

Apr 

13.5  x I02 

4.5  X 10" 

1.34 

.45 

608 

May 

9.4  X 10* 

4.9  X 10* 

2.09 

.60 

S64 

Jun 

4.6  X 102 

4.8  X 10" 

4.17 

.82 

377 

4.69  X 10' 


Hence,  for  a collector  area  of  40.000.  The  price  of  the  solar 

energy  system  compares  favorably  to  the  cost  of  the 
PWFC,  = MFg=(19.72K40.3  X 103)  conventional  system. 


PWFC  = $.744  X 10s 

S 

Step  4 

Determine  B/C,  where, 

PWFC  -PWFC 
B/C  = 5 


For  an  Ap  of  4 X I04  sq  ft 


CWF  = 3.42  X I04  + 10.46  A 


= 3.42  X I04  + 10.46(4  X 104)=S4.53X  105 


_ 1 .38  X 10G  - .794  X 106 
7 " 4.53  X 106 

B/C  = 1 .29 


Table  C4  shows  results  of  the  analysis  for  othei  col- 
lector areas.  A maximum  bene  fit -to-cost  ratio  is  seen  to 
exist.  At  this  point  the  computer  simulation  program 
should  be  consulted  for  a mote  exact  analysis. 


Table  C4 

Benefit-to-Cost  Katin 
for  Various  Collector  Areas 


Collector  Area 
(sq  ft) 

B/C 

0 

0 

5.000 

1.35 

15,000 

1.50 

40,000 

1.29 

85,000 

.99 
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